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The  Air  Force  Office  of  Scientific  Research  (AFOSR)  contract  is  concerned 


with  experimental  and  theoretical  investigations  on  electrically  conducting 


aerosols  for  electromagnetic  obscuration  applications.  Presented  below  is  the 


progress  made  in  the  theoretical  analysis  of  aerosols  consisting  of  magnetic 


(Fe,Co,Ni)  particle  chains  in  air.  Particular  attention  is  given  to  the 


effects  of  an  external  magnetic  field  on  the  chain  equilibrium  and  the 


coagulation  dynamics. 


1.  A  theory  of  the  diffusion  of  macroscopic,  magnetic  particles  (sus¬ 


pended  in  gaseous  or  llqul^  media)  in  density  and  magnetic  field  gradients  has 


been  developed  from  first  principles  (Fokker-Planck  equation),  in  which  the 


Influence  of  the  random,  fluctuating  magnetic  fields,  produced  collectively  by 


the  magnetic  dipole  particles  in  their  thermal  motions  on  the  diffusing  parti¬ 


cle,  is  taken  into  account  in  a  self-consistent  way.  As  an  application,  the 


steady-state  boundary-value  problem  for  the  diffusion  of  magnetic  grains  in 


the  Inhomogeneous  magnetic  field  of  an  adsorbing  sink  dipole  and  an  external, 


homogeneous  magnetic  field  has  been  solved  by  means  of  a  stream  and  Green's 


functions  approach.  The  coagulation  coefficient  for  magnetic  dipole  particles 


in  the  presence  of  an  external  magnetic  field  have  been  derived  and  applied  to 


the  coagulation  of  magnetic  grains  and  the  formulation  of  magnetic  chains  in 


magnetic  aerosols  with  an  external  magnetic  field. 
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2.  Idealized  statistical  calculations  of  chaining  in  a  dilute  suspension 
of  macroscopic  magnetic  particles  in  a  rarefied  gas  have  been  made  in  the 
presence  of  an  external,  homogeneous  magnetic  field.  The  primary  colloidal 
particles  are  assumed  to  be  spherical,  of  the  same  size,  and  have  saturated 
magnetic  moments.  The  magneto-chemical  potentials  and  the  association- 
dissociation  equations  have  been  derived  for  chains  consisting  of  v  >  1  mag¬ 
netic  grains  as  a  function  of  the  temperature  T,  the  density  Nv  of  chains, 
and  the  homogeneous  magnetic  field  B0.  High  field  Intensities  B0  are 
shown  to  shift  the  chain  length  distribution  F  ■  F(v)  in  favor  of  long  chains, 
v  >  1,  whereas  increasing  temperatures  T  move  the  maximum  of  this  statistical 
distribution  to  smaller  chain  lengths,  v  1.  The  theory  appears  to  be  in 
qualitative  agreement  with  experiments  using  an  external  magnetic  field  for 
the  alignment  of  chains  and  their  stabilization. 

i 

3.  The  artificial  charging  of  magnetic  aerosols  is  of  interest  with 
regard  to  the  reduction  of  the  coagulation  rate  of  Coulomb  repulsion  of  the 
magnetic  particles.  For  this  reason,  the  interesting  properties  of  charged 
particle  gases  and  their  fluctuating  electric  microfields  were  studied.  It  is 
shown  that  (1)  the  collective  electric  fields  act  at  distances  larger  than  the 
characteristic  repulsion  distance  D  ■  (KT/4wne2) */2  of  like-charged 
particles  and  (ii)  the  average  collective  electric  field  is  Ey  - 
( 1 2wnKT )  *  /  2  for  ideal  gases  of  particles  of  the  same  charge  e.  Thus,  in  a 
thermal  equilibrium  gas  of  like-charged  particles,  the  longitudinal 
mlcroflelds  are  considerably  stronger  than  in  a  plasma,  since  in  the  latter 
the  random  electric  fields  of  the  negative  and  positve  charges  nearly 


compensate  each  other.  In  addition,  the  interrelation  between  (average)* 
kinetic,  interaction,  collective  field,  and  electric  self-energies  has  been, 
calculated  for  charged  particle  gases. 

4.  Finally,  an  investigation  on  the  reformulation  of  electromagnetic 
theory  for  space  applications  was  carried  through  in  which  the  existence  of 
(i)  a  preferred  frame  of  reference  and  (ii)  an  electromagnetic  carrier  for  the 
electromagnetic  waves  is  taken  explicitly  into • account.  In  view  of  the  exper¬ 
imental  discovery  of  the  cosmic  microwave  background,  (1)  a  preferred  frame  of 
reference  and  (11)  an  electromagnetic  substratum  can  no  longer  be  denied.  The 
latter  work  was  supported  only  in  part  by  the  AFOSR. 

A  list  of  four  papers  follows.  Two  of  them  have  been  accepted  for  publi¬ 
cation,  while  the  others  ye  still  under  review.  Ue  ask  the  AFOSR  to  accept 
reprints  of  these  publications  as  a  final  report.  These  papers  are  attached 
as  Appendices  A  through  D,  respectively. 
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APPENDIX  A 


DISSOCIATION-ASSOCIATION  EQUILIBRIUM  OF  MAGNETIC  PARTICLE  CHAINS 

IN  HOMOGENEOUS  MAGNETIC  FIELDS 

H.  E  Wilhelm 
Michelson  Laboratory 

Naval  Weapons  Center,  China  Lake,  CA  93555 
ABSRACT 

An  idealized  statistical  theory  of  chaining  in  a  dilute  suspension  of  macroscopic  magnetic 
particles  in  a  rarefied  gas  is  presented  when  an  external  homogeneous  magnetic  field  is  present. 
The  primary  colloidal  particles  are  assumed  to  be  spherical,  of  the  same  size,  and  to  have  saturated 
magnetic  moments.  The  magneto-chemical  potentials  and  the  association-dissociation  equations 
are  derived  for  chains  consisting  of  v  a  1  magnetic  grains,  in  dependence  of  the  temperature  T,  the 
density  N„  of  chains,  and  the  homogeneous  magnetic  field  B0.  High  field  intensities  B0  are  shown 
to  shift  the  chain  length  distribution  F  =  F(v)  in  favor  of  long  chains,  v  1,  whereas  increasing 
temperatures  T  move  the  maximum  of  this  statistical  distribution  to  smaller  chain  lengths,  v  -»  1. 
The  theory  appears  to  be  in  qualitative  agreement  with  oven  experiments  using  an  external 
magnetic  field  for  the  alignment  of  the  chains  and  their  stabilization. 


I 


I  INTRODUCTION 

It  is  known  that  external  electric  or  magnetic  fields  shift  the  molecular  dissociation 
equilibrium  in  favor  of  the  paraelectric,  respectively,  paramagnetic  reaction  products.1 
Statistical  considerations  show  that  a  homogeneous  magnetic  field  B0  affects  the  thermal 
ionization  equilibrium  a  *5  i  +  e  of  a  plasma  if  the  quanta  of  the  oscillatory  electron  (e) 
motion  in  the  B0  field  that  corresponds  classically  to  the  electron  gyration  are  of  the  order  or 
larger  than  the  thermal  energy,  Au  S  KT,  where  u  =  eB<>/m  (e  is  the  charge  and  m  is  the 
mass  of  an  electron).2  The  ionization  may  be  enhanced  or  depressed  depending  on  whether 
the  atoms  (a)  are  diamagnetic  or  paramagnetic.2  The  electron  (eJ  and  positron  (e  +  )  densities 
of  the  thermal  vacuum  equilibrium  e+  +  e_  £*  2y  (gamma  quanta)  are  increased  by  a 
homogeneous  magnetic  field  B0  (in  favor  of  the  paramagnetic  or  spin  particles  e  +  )  3 
Similarly,  the  electron-hole  equilibria  in  solids  are  shifted  towards  higher  electron  and  hole 
concentrations  by  a  hor-  >geneous  magnetic  field,  in  particular  in  crystal  structures  with 
small  effective  electron  and  hole  masses.4  However,  these  interesting  physical  effects  are  in 
general  quantitatively  not  very  significant  since  the  energy  pB0  of  the  magnetic  moments  p  in 
the  magnetic  field  is  small  compared  with  the  thermal  energy,  pB0  <*  KT,  except  in  the  case  of 
suprathermal  fields  B0  2r  KT/p,  which  can  be  generated  through  magnetic  flux 
compression.5** 

Magnetic  field  effects  on  reaction  equilibria  are  quantitatively  extremely  important  if 
the  reacting  species  are  macroscopic  or  colloidal  particles  that  have  paramagnetic  or 
ferromagnetic  moments  p  pg,  which  are  large  compared  with  the  Bohr  magneton  pe  = 
eA/2m  =  9.274  X  10'S4  Am*.  (An  analopeous  conclusion  holds  for  paraelectric  or  ferroelectric 
macroparticles  in  electric  fields.)  For  the  latter  reason,  the  magnetic  dipole  energy  pB0  may 
already  be  larger  than  the  thermal  energy,  pB0  >  KT,  for  moderate  magnetic  field  strengths 
B0,  so  that  the  reaction  equilibrium  is  significantly  affected  by  the  magnetic  field.  In 
particular,  this  is  true  for  colloidal  suspensions  of  saturated  paramagnetic  or  ferromagnetic 
grains  of  radius  a  -  lO-*  m,  which  associate  to  long  chains  consisting  of  v  i>  1  grains  in  the 
presence  of  an  external,  homogeneous  magnetic  field.  Dispersons  of  such  colloidal  chains  of 
electrically  conducting,  ferromagnetic  particles  (Fe,  Co,  Ni)  in  the  atmosphere  are  of  technical 
interest  as  wide-band  electromagnetic  obscurants 

Herein,  we  analyze  the  association-dissociation  equilibrium  of  chains  aN  consisting  of 
v  spherical  grains  of  radius  r  =  a  and  saturated  magnetic  moment  pi  dispersed  in  dilute. 


homogeneous  carrier  gases  of  temperature  T  and  a  homogeneous  external  magnetic  field  B0.  • 
The  corresponding  reaction  equations  are 

a  t*  a  ,  +  a,  ,  2  £  v  <  v  ,  d) 

where  &i  *  a  and  v»  <  ®  is  the  maximum  number  of  grains  in  a  chain.  If  the  grain  material 
is  treated  as  a  single,  magnetic  domain  of  magnetization  M[A/m),  the  magnetic  moment 
[  Am2)  of  a  spherical  grain  of  radius  r  =  a  is  given  by 

p  =  4nm/6ki3M  .  (2* 

In  local  magnetic  field  B0  the  potential  energy  of  a  grain  is  U  =  -p  •  B0  ^  0.  The  average 
moment  <p  | )  of  a  grain  in  direction  of  field  B0  is  in  thermal  equilibrium  (p  | )  =  L(pB<j/KT) 
where  Lie,,)  =  cothe0  -  e„'‘  =*  1,  e0  >  1,  is  the  Langevin  function.  Accordingly,  <p  | )  =  p  for 
pB0  i>  KT.  The  maximum  potential  energy  of  the  dipoles  in  the  magnetic  field  B0  (relative  to 
the  thermal  energy'  KT)  is  by  Eq.  (2) 

c  =  4n(n/6)a3MB/KT  =  4.77  x  1023 a3MBlT  <3> 

O 

Accordingly,  c0  -  1 0*  for  M  =  10*/4n  A/m,  B0  =  1  Vsec/m2,  T  =  300°K,  and  a  =  lO-”  m.  This 
example  shows  that  at  standard  temperatures  T,  the  grains  have  their  moment  p  aligned 
parallel  to  B0  already  at  moderate  field  intensities,  (p_'  =  p  H  B0  for  e  *>  1. 

Since  the  magnetic  field  of  a  dipole  pt  is  Bi  =  -V^pj  •  r/4nr’l,  the  interaction  energy 
of  two  magnetic  dipoles  ”i"  amd  "j"  at  a  distance  r  =  j  r  |  apart  is  Uy  =  (po/4n)r'3[p,  •  pj  - 
Sr^lp,  •  r)(pj  •  r)l  where  m,  =  4n  X  10'T  Vsec/Am.  The  average  distance  between  grains  is  on 
the  order?  -  n  i3  where  n(m  31  designates  their  density  Accordingly,  the  average  binary 
interaction  energy  of  the  dipoles  is 

c  =  4np  (n/6)2a6M2n/KT  =  3.14  X  1017a6M2n/T  (4) 

IJ  r0 

relative  to  the  thermal  energy  KT.  Hence,  T,j  =  6.63  x  10'’  <5  1  for  M  =  10f/4n  ,Vm,  T  = 
300°K,  a  =  10"*  m,  and  n  S  10*'  nr*.  Thus,  the  grains  behave  like  an  ideal  gas  atT  =  300°K 
for  densities  n  <  1021  m  J  In  the  following  it  will  be  assumed  that  the  grains  and  the  chains 
formed  from  them  behave  like  a  quasi-ideal  suspension  in  the  atomic  (nonmagnetic)  carrier 


gas.  By  Eq.  (4)  this  requires "c, j  <<  1,  i  e.,  sufficiently  low  particle  densities  n  and  sufficiently  • 
high  system  temperatures  T  (below  the  Curie  point). 

Physically  more  complex  is  the  calculation  of  the  chain  lengths  equilibrium  for 
ferromagnetic  grains  dispersed  at  highest  densities  (n  -  1023  m'3)  in  so-called  ferromagnetic 
fluids.1  For  this  situation,  an  asymptotic  theory  for  strongly  nonideal  interactions  was 
proposed  (which  breaks  down  in  the  ideal  limit).7  This  nonideal  theory  has  been  shown  to 
disagree  up  to  orders  of  magnitude  with  the  experimental  data.8 


*.J.  4-> 


y  M.  a-t. 


II.  MAGNETOACTIVE  CHAIN  EQUILIBRIUM 


A  homogeneous  magnetic  field  B0  aligns  the  magnetic  dipole  moments  pi  of  the 
grains  at  and  thus  provides  favorable  conditions  for  the  formation  of  linear  particle  chains 
resembling  rather  rigid  one-dimensional  polymers.  In  thermal  equilibrium,  the  chain 
lengths  Ly  or  the  numbers  v  =  1,  2,  3,  ....  v»  of  grains  at  in  the  v-chains  have  a  distribution 
F  =  F(v)  determined  by  the  maximum  entropy  principle.9  The  most  elementary  dissociation 
(-*)  -association  («-)  reactions  are  those  in  which  one  grain  at  is  removed  or  added  to  a 


v-chain  from  the  collective  2  S  v  <  v,: 


a  *5  a  ,  +  a.  , 
v  v-t  l 


a  * 5  a  „  +  a,  , 
v  v-2  1  ’ 


°3t5a2  +  fll  ’ 


a2*5aI  +  a| 


Adding  these  relations  yields  a  summary  reactin  equation,  which  describes  the  dissociation  of 
a  v-chain  into  v  grains  at  and  the  formation  of  a  v-chain  by  association  of  v  grains  at, 
respectively: 


\-val  , V  5=  2, 3, 4 - Va 


This  reaction  is  of  particular  interest  since  it  relates  the  densities  Nv  of  the  chains  v  >  2  to 
the  density  N  |  of  grains  at .  The  distribution  Nv  of  the  v-chains  observed  in  experiments  is  the 
one  with  largest  probability  and,  hence,  maximizes  the  entropy  of  the  colloidal  suspension 
The  latter  condition  leads  to  a  statistical  equilibrium  equation  for  the  summary  reaction  (6) 
in  terms  of  the  magneto-chemical  potential  Qv  of  the  v-chains  and  the  magneto-chemical 
potential  Q]  of  the  grains  aj. 


Qv  =  vQj  ,  v  =  2, 3, 4,  . .,  ■ 


The  magneto-chemical  potentials  of  the  chains  v  >  1  and  the  grains  v  =  1  are  given  for  ideal 


conditions  bv 


•>» 


A  A'A*  V  »  V 


where 


- Q  =  KTilnlN^)  +7JKT  -  ln(U™U8)} 


(8)  • 


A  =  hJMnm  KDl/2  ,  (9) 

V  V  9 

UVR  =  u v  UR  ,  ,10) 

V  V  V 

UB  =  sinh(p  B  IKT)/(p  B /KT)  ,  <”> 

t  V  r  V  o  r  v  o 

I 

i  are  the  thermal  OeBroglie  wavelength  of  the  v-chains  of  mass  mv,  the  product  of  their 

vibrational  (V)  and  rotational  CR)  partition  functions,  and  the  orientational  partition 
functions  of  the  magnetic  moments  pv  of  the  v-chains.  Obviously, 


J2n  in 

d<J>  j  expip cos  8)  sin  9  d8  , 


(12) 


where  <J>  =  <  (pv>  B0).  The  masses  mv  and  moments  pv  are  given  in  terms  of  mi  and  pi, 
respectively, 


n  =  vm.  .  pu  =  vp 


(13) 


Substitution  of  Eq.  (8)  into  Eq.  (7)  yields  the  fundamental  dissociation-association 
equations  for  v-chains  in  a  homogeneous  magnetic  field  B0: 


p,  sinh(p  B  / KT)  i  pB  u 
JV  o  l"lo  * 

v  1  '  v 


PvIsin/ifpjB^/fDr 


KT 


(2nml/f7')' 


,3(2 


v- 1  /Ac 

exp' 


(*) 

'KT/ 


v  =  2,3,4.  ,,v 


114) 


where 


A?  ■  vc,  -T  ~  v(poPj/4n(2o1)3)  , 


fVR  m  uVRnuVRii  -  10° 

’  V  V  I 


(15) 

(16) 


The  energy  Aev  required  for  the  dissociation  of  a  v-chain  into  v  grains  is  of  the  magnitude 
Ac^  -  vlO'1  eV.  The  vaJues  of  A?v  can  be  calculated  (v-body  problem)  or  measured 
experimentally.  Equation  (16)  holds  in  view  of  the  macroscopic  nature  of  the  particles  v  ^  1 

Equation  (14)  represents  v»  -  1  equations  for  the  determination  of  the  v*  unknown 
particle  densities  Nj,  NT2,  N3, .  .  NTVv  A  complete  system  of  equations  is  obtained  by  adding 
the  equation  for  the  conservation  of  the  total  grain  density  N, 


y  vAf  =  N 

V 

V  *  1 

Substitution  of  Eq  (14)  shows  that  Eq.  (1?)  represents  a  polynomial  of  order  vx  for  the 
determination  of  the  grain  density  N’i, 


V  C(B  =  N 

V  o'  1 


where 


p  sinh(p  B  IKT)  p.B  f,3  v-i  /Ac  • 

C  (B  ,T)  =  vT“ - — -  - - -  exp(  —  , 

v  °  v  PAsinh(pxBJKT)\ v  * T  (inmfiTf1  '  KT  ' 

v  =  2,3,4 . 


With  N't  calculated  as  the  positive  root  of  Eq.  (18),  the  chain  densities  follow  from  Eq  (14)  as 
N’v  =  CV(B0,  T)Nj  v 

In  applications,  one  is  mainly  concerned  with  suprathermal  magnetic  fields  for  which 
B  i>  KT/p ,  :  p  B  /KT  P  1  .  (20) 

O  r  1  *V  0 

For  such  strong  magnetic  fields,  Eq  (19)  reduces  to  the  dissociation-association  equation 


AT  =  N'i2s‘~\V2fVR 
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v-i  /Ac 

«p 


(“I 

'  KT  )  ’ 


v=  2,3.4 _ v  S>#C7Vn,  , 

*  u  M  ’ 
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KWWWW 


with 


C  (B  1D  =  2v',v1/2/,V/? 

VO  9  V 


P,B 
^  1  o 


KT  ttnm.KT)3'2 


v- 1  /  Ac 

exp 


/Ac  x 

(«?)• 


=  2,3,4,..  ,v.  ,Boi>KT/pl  . 


(22) 


Equation  (21)  provides  interesting  proportional  relations  that  exhibit  the  main 
magnetic  field  (B0)  and  temperature  (T)  dependences  of  the  chain  densities  Nv: 


N2«  N2]B1T~s/2  exp[\tJKT) 


N3  «  Mffir  ~ 1<y2  exp^yKT) 


Na  «  AT JflJr  ' 1W2  expdSjKT)  ,  Bg  *>  KT/p, 


(23) 


NV«  AfjB*"1T*5,v‘lW2  exp(teJKT) 


It  is  seen  that  an  increase  of  the  magnetic  field  B0  shifts  the  distribution  {Nj,  N2, . . NVJ  in 
favor  of  the  large  chains,  v  -*  v*,  whereas  an  increase  in  the  temperature  T  shifts  the 
distribution  {Nv}  in  favor  of  the  short  chains,  v  -*  l.' in  view  of  the  exponential  T-dependence 
in  Eq.  (23)  and  Acv  >  Acv_i,  a  temperature  increase  has  a  particularly  strong  destructive 
effect  on  the  long  chains. 


V-V 


The  statistical  theory  presented  should  be  considered  a  first  step  towards  the 
qualitative  understanding  of  the  dissociation-association  equilibrium  of  magnetic  colloid 
chains  in  an  external,  homogeneous  magnetic  field,  presumed  that  the  primary  grains  formed 
originally  a  sufficiently  dilute  suspension  in  a  nonmagnetic  gas.  The  favorable  effect  of  the 
external  magnetic  field  on  the  chain  formation  and  the  shifting  of  the  chain  lengths 
distribution  F(v)  towards  larger  chain  lengths  v  has  been  observed  experimentally.  However, 
quantitative  experimental  data  on  magnetic  chain  lengths  distributions  in  aerosols,  which 
could  be  used  for  comparison  with  the  theory,  have  apparently  not  yet  been  published.  The 
experiments  also  indicate  that  chains  and  grains  coagulate  to  large  clusters  which  "fall  out” 
in  the  gravitational  field  when  the  magnetic  field  is  switched  off  The  rapid  coagulation  in  the 
absense  of  a  (sufficiently  strong)  magnetic  field  can  be  reduced  by  spraying  the  macroscopic 
particles  with  electric  charges  (Coulomb  repulsion).*® 

In  order-  to  achieve  a  quantitative  understanding  of  the  chain  equilibrium  in  an 
external  magnetic  field,  several  (difficult)  problems  remain  to  be  solved  rigorously.  The 
dissociation  energies  Ac"v  =  vci  -Tv  are  to  be  calculated  from  the  v-body  problem  of  a  v-chain 
made  up  of  extended  magnetic  dipoles  (grains  of  radius  a  >  0).  Based  on  this  v-body 
dynamics,  the  oscillatory  and  rotational  partition  functions  of  the  v-chains  have  to  be 
determined  (even  though  these  degrees  of  freedom  of  macroscopic  particles  are  only  poorly- 
excited  at  temperatures  below  the  Curie  point).  Finally,  to  render  the  theory  applicable  to 
higher  grain  densities,  the  nonideal  interactions  between  chains  and  grains  would  have  to  be 
taken  into  consideration  in  the  evaluation  of  the  statistical  chain  equilibrium. 

In  connection  with  the  chaining  phenomenon,  various  other  interesting  effects  could 
be  investigated  theoretically.  For  example,  it  would  be  important  to  understand  the  effect  of 
primary  grains  that  have  not  the  same  size  but  a  size  distribution.  Furthermore,  it  would  be 
interesting  to  evaluate  the  chain  lengths  distribution  for  primary  grains  that  carry  an 
artificial  electric  charge.  Already  these  few  examples  indicate  that  the  research  on 
magnetically  active  colloids  offer  significant  opportunities  for  further  contributions  *® 
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APPENDIX  A:  Alternative  Derivation 


It  is  instructive  to  derive  the  dissociation-association  equation  (Eq.  (14)1  also  from  the 


reaction 


where 


a  **  a  ,  +  a,  , 

v  v-l  1 

v  =2,3,4... 

>v»  • 

(A-l) 

fi  =  Q  .  +  0.  . 

V  v-l  1 

v  =2,3,4... 

-v.  - 

(A-2) 

in  statistical  equilibrium.  Substitution  of  the  magneto-chemical  potentials  [Eq.  (8)1  into  Eq. 
(A-2)  yields  (after  some  alegebra)  the  "one-grain”  dissociation-association  equation: 


where 


IV  =D  (B  ,  T)N  N  ,  v  =  2,3,4 . v  , 

V  Vo  V-l  1  9  *  00  » 

DIB  n_W  v  )"*- 

v  °'  v  '  v  —  1  /  p  sin ht p  ,B  !KT)sinh 


sinhlpv_  XBJK T)si nh(p} BjK T) 


plBo  h3  _  _  1 

x  ™  exp  (T  +e  -c  V/iCT  ,  v  =  2,3,4... ,v 

KT  (2nm1KTTt  11  » 


*VR  =  fjVRHJVR  JJVR  (A-5 

®V  V  v-l  l 

Eq.  (A-3)  represents  a  recurrent  relation  that  gives,  by  elimination,  the  "v-grain” 
dissociation-association  equation: 


where 


•  v  =  2,3,4...,v<b  , 

v=  2 


n  0 ,Bo. r,  =  G„,Knr—  \ ’"  **(  T? ) 


-l  /vz.-T  v  (A-7) 
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and 


since 
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(t  .  +  c 
v- 1  1 


-£v)  = 


VC,- 


( A-9] 


Factorization  of  the  finite  product  (A-8)  and  evaluation  of  the  individual  products  gives 


since 


G  -f  xvwx 

V  V 


Pj  sinhlp  B  IKT) 

—  x  - — - 

pv  sinh{p  B  IKT) 


n«v 


v-2 


(A- 10) 


(All) 


by  Eqs.  (A-5)  and  (16).  Combining  Eqs.  (A-6),  (A-7),  and  (A-9)  results  in  the  dissociation- 
association  equation: 


N 


V 


jv;Ak 

l' V 


Pt  sinhipJBjKT) 
Pvlsinh{p{Bo/KT)}v 


M.  »3  I-1  /vV7,\ 

KT  aunjerP  I  KT  I 


V 


2,3,4 


(A-12) 


Equation  (A-12)  is  identical  to  Eq.  (14)  derived  from  the  summary  reaction  [Eq.  (6)1  and  the 
statistical  equilibrium  relation  (7) .  The  derivation  of  Eq.  (14)  is  apparently  simpler  than 
that  of  Eq.  (A-12).  The  identity  is  necessary  for  physical  reasons. 
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ABSTRACT 

A  theory  of  the  diffusion  of  macroscopic,  magnetic  particles  (suspended 
in  gaseous  or  liquid  media)  in  density  and  magnetic  field  gradients  is 
developed  from  first  principles  (Fokker-Planck  equation).  The  influence  of 
the  random,  fluctuating  magnetic  fields,  produced  collectively  by  the  magnetic 
dipole  particles  in  their  thermal  motions  on  the  diffusing  particle,  is  taken 
into  account  in  a  self-consistent  way.  It  is  shown  that  the  anisotropy  in  the 
particle  diffusion,  caused  by  the  coupling  of  translational  and  rotational 
degrees  of  freedom  (Magnus  effect),  is  small  in  most  physical  situations.  As 
an  application,  the  steady-state  boundary-value  problem  for  the  diffusion  of 
magnetic  grains  in  the  inhomogeneous  magnetic  field  of  an  adsorbing  sink 
dipole  and  an  external,  homogeneous  magnetic  field  is  solved  by  means  of  a 
stream  and  Green's  functions  approach.  The  coagulation  coefficient  for 

magnetic  dipole  particles  in  the  presence  of  an  external  magnetic  field  is 

* 

derived.  The  results  are  discussed  with  regard  to  the  coagulation  of 
magnetic  grains  and  the  formation  of  magnetic  chains. 


INTRODUCTION 


The  physical  behavior  of  colloids  formed  from  ferromagnetic  grains 
suspended  in  magnetically  Inactive  liquids  has  recently  been  investigated 
experimentally1*2  and  theoretically.3*4  Under  the  influence  of  a  uniform 

external  magnetic  field,  the  magnetic  grains  tend  to  associate  in  the  form  of 
long  chains  if  the  grain  density  and  the  temperature  of  the  liquid  are 
sufficiently  low.5*6 

Based  on  the  Fokker-Planck  equation,  we  develop  an  analytical  theory  of 
the  thermal  diffusion  of  magnetic  dipole  particles  in  magnetically  passive 
fluids  F  (gases  or  liquids)  when  an  inhomogeneous  magnetic  field  TfdT.t)  is 
present.  A  novel  diffusion  equation  is  derived  and  applied  to  the  diffusion 
of  magnetic  grains  in  the  inhomogeneous  field  B$(r)  of  an  ’’absorbing" 
magnetic  dipole  when  a  homogeneous  external  field  is  present.  The 
corresponding  boundary-value  problem  is  solved  analytically.  The  coagulation 
of  dipole  particles,  which  interact  through  their  magnetic  self-field,  is 
discussed  in  dependence  of  the  external  field  B^. 

Ferromagnetic  grains  have  a  typical  radius  a  —  10~8  m.  If  the  grain 

material  la  treated  as  a  single  magnetic  domain  of  magnetization  M[ A/m] ,  the 
magnetic  moment  [Am2]  of  a  spherical  grain  of  radius  r  ■  a  is  given  by 

p  ■  4ir(ir/6)a3i?  .  (1) 

Let  these  grains  be  dispersed  in  a  nonmagnetic  (p0)  carrier  medium  (F)  of 

temperature  T.  In  a  local  magnetic  field  T,  the  potential  energy  of  a  grain 

1®  +  *  s 

U  -  -  p*B  $  0  (2) 


The  average  moment  <pn>  of  a  grain  in  direction  of  field  B  is  in  thermal 
equilibrium 


<P|>  -  L(e)p,  e  =  pB/KT  ,  (3) 

where 

L(e)  "  cothe  -  e"1  ■  1,  e  >  1  (4) 

is  the  Langevin  function.  Accordingly,  <pu>  *  p  for  pB  >  KT.  The  maximum 
potential  energy  of  the  dipoles  in  field  T  (relative  to  the  thermal  energy  KT) 
is  by  Eq.  (1) 

e  -  4n(w/6)a3MB/KT  -  4.77  x  1023  a3MB/T  .  (5) 

Accordingly,  e  —  102  for  M  ■  10®/4ir  A/m,  B  *  1  Vsec/m2,  T  ■  300°K,  and  a  ■ 
10“8  m.  This  example  shows  that  at  standard  temperatures  T,  the  grains  have 
their  moment  ~p  aligned  parallel  to  B  already  at  moderate  intensities,  <jT«>  * 
~p  *  if  for  c  >  1. 

Since  the  magnetic  field  of  a  dipole  "pi  is  ?i  “  -P(uoPi  *T/4irr3) , 
the  Interaction  energy  of  two  dipoles  "1"  and  "j”  at  a  distance  r  •  |7j  apart 
Is  (u0  *  4ir  x  10”7  Vsec/Am) 

Uij  -  (u0/4ff)r"3tPl*Pj  “  3r"2(Pi»r)(Pj .?)]  .  (6) 

The  average  distance  between  the  dipoles  is  of  the  order  r  —  n-1^3  where 

n[m~3]  designates  their  density.  The  average  binary  Interaction  energy  of  the 
dipoles  is  by  Eq.  (6)  of  the  order 

eiJ  “  4iUi0(»/6)2a6M2n/KT  "  3*14  *  i017a6M2n/T  (7) 

relative  to  the  thermal  energy  KT.  Hence,  ej.j  *6.63  x  10" 3  <  1  for  M  - 

106/4»  A/a,  T  •  300°K,  a  •  10"  a  a,  and  n  <  1021  m"3.  Thus,  the  grains  behave 


-2- 


teJ&j*** 


■  v-vywi.- >•  v- . 


*  mVlC%**_**P 1 


t  *>  V  ’  •  ‘ 


like  an  Ideal  gas  at  T  *  300°K  for  densities  n  £  10 21  m“3.  The  raaximun  binar;; 
interaction  energy  of  the  dipoles  is  (r  ■  2a) 

“  u0*(*/6)2M2a3/KT  -  7.84  x  1016M2a3/T  (8) 

by  Eq.  (6).  Accordingly,  ejj  —  10°  for  M  *  106/4it  A/m,  T  *  300°K,  and  a  * 
10“8  m.  Therefore,  under  typical  conditions,  some  of  the  grains  will  always 
coagulate  as  a  result  of  their  dipole  attraction. 

These  numerical  illustrations  Indicate  that  the  grains  coagulate  in  their 
attracting  inhomogeneous  dipole  fields  Bs  to  larger  macroparticles,  in 
particular  at  low  temperatures.  In  homogeneous  external  fields  of  moderate 
Intensity,  B0  <  1  Vsec/m2,  the  coagulation  leads  to  long  chains  consisting 
of  up  to  106  grains  which  are  aligned  parallel  to  B0.  In  thermal 
equilibrium,  the  chain  length  distribution  is  determined  by  the  respective 
temperature  T  and  grain  density  n. 


ta  Jm  >.•  fa'  a  J  i  a*  1*  *at,fc>,  ta^  A*  A*  V  M  Jl‘ J^*  Jl*  >i*  hf  !■*  Ih»  fa*  N*  fa*  >«■»  1w*  fa*  ».1.L*  VI  »,*  u  U  «.♦  tm'tjf  K*  «*■ 


KINETIC  EQUATION 

The  diffusion  of  macroscopic  ferromagnetic  grains  of  radius  a,  mass  -m, 
moment  of  inertia  I,  angular  momentum  L,  and  magnetic  moment  7  in  a 
nonmagnetic  carrier  medium  F  (gas  or  liquid),  and  an  inhomogeneous  magnetic 
field  If  represents,  in  general,  a  physically  complex  problem.  The  main 
complications  are  (i)  the  coupling  of  translational  (v)  and  angular  ( .j) 
velocities  of  each  grain  (Magnus  effect)  and  (ii)  the  precession  of  the 
magnetic  moment  about  the  locaL  magnetic  field  3  in  case  of  incomplete 
alignment.  The  dynamical  variables  L,  oj ,  and  p  of  a  grain  are  interrelated  by 


L  »  loi. 


where 


df/dt  •  p  x  J,  dp/dt  *  yp  *  $  (10) 

determine  the  precession  of  L  and  "p  in  8.  The  gyromagnetic  factor  is  y  =• 
p0g/f»  (8  “  Landd  ratio,  pQ  *  efi/2me  »  Bohr  magnetron),  since  the  ferro¬ 
magnetism  of  the  grains  is  due  to  their  electrons.  The  precession  of  L  and 
about  the  direction  of  B  occurs  with  Che  modified  Larmor  frequency7 

•  -wfl,  ■  egB/Zme  >  0  (11) 

In  view  of  Eq.  (9),  L,  p,  and  u  are  equivalent  dynamical 1  variables  for 
any  grain.  Th«  distribution  function  of  the  translational  ~v  and  angular 
velocities  la  of  the  grains  is,  therefore,  at  any  point  (r”,t)  of  the  form 


f  ■  f (r, t,v,u) 


where 


are  the  macroscopic  moments  (grain  density  n,  grain  flux  n<v>,  ...)  of 

f  ( r ,  t  .  Changes  of  the  distribution  function  are  brought  about  by  the 

force  of  the  inhomogeneous  B  field  on  the  magnetic  moment  ~p,  the  translational 


friction  force  (relaxation 

time  Tj), 

the  Magnus  force, 

and  the  friction  torque 

(relaxation  time  r2)  on  the 

grains: 

^m  "  P 

(14) 

HI 

fi 

1 

mv/ tj  , 

(15) 

?c  - 

icmav  x  a)  , 

(16) 

«f  ■- 

it/ t2  . 

(17) 

The  equivalent  fluid  mass  of  the  grain  volume  is  ma  *  ^na^mp/B. 
The  Magnus  force  (Eq.  (16)]  couples  the  translational  ~v  and  angular  motions 
of  each  grain  such  that  a  grain  rolls  aside  if  T  is  not  parallel  to  oT.  8»9 
The  dimensionless  correlation  Integral  of  ~:c  is  of  the  order9  tc(f()  ~  10°. 
The  frictional  relaxation  frequencies  of  T  and  "w  are  in  the  free  molecular 
region,10  K  -  X^/a  >1: 

T”1  -  (8/3)a2np(2innFKT)1,2(l  +  aiff/8)/m  ,  (18) 

T"1  -  (32/9)a4nF(27rmFKT)l/2(l  +  o2it/8)/I  , 


(19) 


where  0  <  a<  j  <  1  are  accommodation  coeff icients 10  (m_  ■  mass,  n  = 

density  of  carrier  fluid).  In  gases  with  a  m.f.p.  Xp  ~  10~7  m,  the 

Knudsen  number  is  K  ~  10 1  for  grains  with  radii  a  ~  10~8  m.  The  corresponding 

Stokes  formulae10  (n  “  viscosity  of  carrier  fluid)  hold  in  the  continuum 

F 


region,  K  -  Xp/a  «  I : 


r" 1  *  6irn_a/ra,  r  1  =*  Sirr^a3/! 


For  obvious  physical  reasons,  the  distribution  function  f(£*, t,~v7uj)  of  the 
grains  satisfies  a  continuity  equation  in  the  {t.T.V.Hr}  space.  Consideration 
of  the  forces  (14)-(16)  and  the  torque  (17)  on  the  balance  of  grains  which 
leave  and  enter  the  volume  element  dV  •  d3r”  d3v"  d  3£T  at  the  point  in 

the  time  dt  ♦  0,  yields  the  fundamental  kinetic  equation  for  the  distribution 
function  f  ( t , v’,~u) : 

3f  3  ♦  3  p»7^  3  ♦ 

—  ♦  —  •  (vf)  +  —  •  ’  f)  +  T  •  (Sf)  - 

3t  3r  3v  3ui 


K  J"(vxu)  •  ♦  +  T”1  ♦  •  [(V  +  m  *)f)  +  T,  ♦  •  Ku>  +  X  *)f]  (21) 

3w  1  3v  3v  2  3u»  3oi 


where 
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7  K  H  -  0 
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72*  -  7.  f  /  p  f  d3v  d3u»  . 


The  Fokker-Planck  equation  (21)  and  the  Poisson  equation  (24)  represent. a 


system  of  integro-dif ferential  equations 

for 

the 

determination 

of 

the 

distribution  function  f(  r”,t  ,7,"^)  and 

the 

potential 

$(T,  t) 

of 

the 

self-consistent  magnetic  field  H(r\t), 

which 

has 

its 

sources 

in 

the 

distributed  dipoles  The  angular  acceleration  is  defined  by  Eq.  (10). 

Equations  (10),  (22),  and  (23)  are  auxiliary  equations,  assuming  the  absence 
of  electric  currents  (V  x  H  ■  0). 

In  Eq.  (21)  the  interaction  terms  with  coefficients  k,  t"1,  and  t2_1  are 
responsible  for  the  relaxation  of  f.  In  an  external  magnetic  field  fy,,  the 
thermal  equilibrium  distribution  is  given  by  (p“ ■  yloD 

f0  "  C0(n0,  To.  Bo)  exp[-(J  mv2  +  j  Iu2  -  p-B0)/KTl  ,  (25) 

If  the  nonideal  -  field  effects  due  to  the  dipole  interactions  are  neglected 
(C0  -  normalization  constant).  By  Eq.  (25),  the  average  energies  of 
translation  and  rotation  are  in  thermal  equilibrium  T0 

<mv2/2>  -  3KT0/2 ,  <Iu2/2>  -  3KT0/2  .  (26) 

For  vanishing  interactions  of  the  grains  with  the  carrier  medium, 
t1(2  *  **,  the  relaxation  terms  in  Eq.  (21)  are  absent.  In  this  case, 
Eq.  (21)  has  elementary,  stationary  (3/3t  s  0)  solutions  of  the  form 


f9  "  H(e), 


1  ♦-  1  ♦  ♦ 
e  =  2  07  +  2  *<*>  ~  /♦(p*VB)  •  dr 


(27) 


3fs/3r  *  -  H'(e)p*VB,  3fs/3v  =  H'(£)mv  , 


3*(wfs)/3uj  »  H '  (  e ) I ( tu  -  •  <j  =■  0  , 


where  u  *  yu  x  B  Is  perpendicular  to  (w  -  yB)  by  Eq.  (10).  In  Eq.  (27),  H(e) 


is  an  arbitrary  functional  of  e  and  H'(e)  =  dH(e)/de.  The  Magnus  term  in 


Eq.  (21)  vanishes  since  3fs/3v  II  v. 


In  general,  nonequilibrium  solutions  of  Eq.  (21)  are  obtainable  by  means 


of  perturbation  theory.  The  latter  assumes  f  *  f0  +  f  for  0  <  t\/2  <  ®, 


with  If  I  €  f0.  It  should  be  noted  that  the  kinetic  equation  (21)  implies 


that  the  grain  component  behaves  like  an  ideal  gas  with  an  equation  of  state 


p  -  nKT. 


DIFFUSION  EQUATION 


In  Eq.  (21),  the  Magnus  force  terra  causes  grains  with  ~7 J-  U  to  diffuse  * 
somewhat  slower  than  those  with  V  )  7.  As  a  result,  the  translational  (v)  and 
rotational  (u)  velocities  are  coupled,  and  the  diffusion  of  grains  becomes 
anisotropic.  As  known  from  the  theory  of  the  Senftleben  effect11"14  for 
molecules  (with  mechanical  spin  and  magnetic  moment)  in  a  homogeneous  magnetic 
field,  the  anisotropy  of  the.  diffusion  coefficent  decreases  with  increasing 
particle  size. 

Nondimenslonallzation  of  Eq.  (21)  Indicates  that  the  strength  of  the 
coupling  of  the  translational  (v)  and  rotational  (w)  motions  by  the  Magnus 
force  «  7  x  J is  determined  by  the  relaxation  frequency 

tc1  -  <(2/3)l/2(ma/m)(3KT/I)l/2  ,  (30) 

where  wp  ■  (3KT/I)1/2  is  the  thermal  frequency  of  rotation,  and  I  »  3raa2/5 
for  spherical  grains.  Accordingly,  one  expects  the  diffusion  anisotropies  to 
be  Insignificant  for  coupling  frequencies  Xj.”1  which  are  small  in  comparison 
with  the  translational  relaxation  frequency  Xj"1. 

The  diffusion  equation  for  the  grains  'follows  from  their  continuity  and 
motion  equations  in  the  carrier  medium.  The  moments  of  the  kinetic  Eq.  (21) 
for  the  corresponding  dynamical  variables  (mv)°  and  (mv) 1  are 


3n/3t  +  V»(n  <v>)  •  0 


(31) 


3(naKv>)/3t  +  V»(nm<vXv>)  ■  -7(nKT)  +  n<p>»7B  -  xj"1  nm(<v>  -  <vp>)  (32) 


for 


Xc1  <  Xl_1(l  +  X2**  V  Xj”  1 )  ^2  xf1,  X2_1  <  Xj 


-1 


(33) 
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Equation  (33)  is  the  condition  for  the  quasi-isotropic  diffusion, 
approximation. 14  The  sources  for  the  momentum  changes  of  the  grains  are  the 
pressure  gradient,  the  magnetic  dipole  force  density,  and  the  intercomponent 
friction  force  density  between  the  grains  and  the  carrier  medium  (F)  with  mean 
mass  velocities  <v>  and  <v^>,  respectively. 

The  magnetic  dipole  force  density  in  Eq.  (32)  results  from  Eq.  (21),  in 
accordance  with  the  partial  integration 


n<p>  »VB 


(34) 


for  particles  with  magnetic  moment  "p  *  yI<T,  and  f  —  exp(-mv2/KT)  ♦  0  for 

For  diffusive  or  slow  shock-free  motions  of  the  grains,  the  nonlinear 
inertia  term  V»(nm<v><v>)  in  Eq.  (32)  is  neliglble.  In  this  approximation, 
elimination  of  n<v>  from  Eqs.  (31)  and  (32)  yields  the  "hyperbolic"  diffusion 
equation  for  the  grain  density  field  n(?,t)  in  isothermal  carrier  media  (T): 

32n/3t2  +  tT1  an/ at  +  t7lV.(n<v  >)-«  c2V» l Vn  -  n(<p>.vS)/KT]  ,  (35) 

*  1  r"  T 

where 

cT  -  (KT/m)1/2  (36) 

is  the  thermal  speed  of  the  grains  of  mass  m.  Equation  (35)  represents  a  wave 
equation  which  propagates  density  perturbations  n(?,t)  with  a  characteristic 
speed  cy  and  relaxation  time  t  j . 

In  applications  usually  only  large  observation  times  are  of  interest, 
t  >  ij.  In  this  case,  Eq.  (35)  can  be  reduced  to  the  "parabolic”  diffusion 


equation  for  grains  in  isothermal  carrier  media  (T): 


3n/3t  +  7»(n<Vp>)  «  D7*f7n  -  n(<p> *7B)/KT] 

where 

Q  *  Tlcx2  *  (KT/mJt! 


(3/;' 


(38) 


Is  the  diffusion  coefficient  In  the  quasi-isotropic  approximation  (33). 

Equation  (37)  propagates  density  perturbations  with  infinite  maximum  speed, 

since  c-j2  ♦  •  for  ♦  0  by  Eq.  (38). 15 

In  Eqs.  (35)  and  (37),  7*(n<v_>)  considers  the  convection  of  grains  by 

the  flow  <v  >  of  the  carrier  medium  (F).  The  self-consistent  magnetic 
F 


field  B  in  these  hyperbolic 

and  parabolic 

diffusion 

equations 

is  described  by 

-► 

♦  +■ 

B  - 

Wo11  +  u0«<P>. 

7*B  -  0 

t 

(39) 

■¥ 

♦ 

♦ 

H  - 

-  7S, 

7  x  H  - 

0, 

(40) 

72* 

»  7«(n<p», 

(41) 

where 

♦ 

♦  ♦ 

d<p>/dt  ■  y<p>  *  B 


(42) 


Equations  (39)-(42)  are  the  macroscopic  (average)  versions  of  Eqs.  (22)-(24) 
and  Eq.  (10),  respectively. 

The  presented  diffusion  theory  for  magnetic  grains  suspended  in  gaseous 
or  liquid  media  (F)  Is  applicable  in  the  isotropic  diffusion  approximation 
(33),  which  requires  that 

Tc1/tJ1  ~  (mF/m) 1/2  «  1,  K  >  1,  (43) 

TcVfJ1  -  (nFmF/nF)(KT/m)1/2a  Cl,  K  «  1, 


(44) 
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k'- 


in  the  free  molecular  and  continuum  flow  regions,  respectively.  For'  K  »  l, 
the  isotropy  condition  is  satisfied  since  €  m  for  macroscopic  grains 
(a  >  10“10  m  )  of  mass  m  *  (4ira3/3)p  (solid-state  density  p  >  103  kg/m3).  For 
K  €  1,  the  istropy  condition  is  satisfied  for  dense  gases,  and  also  for 
typical  liquids  (nFmF  -  iO3  kg/m3,  np  -  10~3  kg/msec,  T  -  300°K)  for 
which  —  10“6/a1/2  «  1  for  a  >  10“ 10  m. 
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BOUNDARY-VALUE  PROBLEM 


In  Che  classical  theory  of  Smoluchowski , 16  noncharged  and  nonmagnetic 
colloidal  particles  of  radii  a  and  b  ^  a  are  shown  to  coagulate  upon 

approaching  their  critical  interaction  sphere  of  radius  d  *  a  +  b  by  Brownian 
motion.  When  the  grains  are  surrounded  by  a  layer  of  solvent,  the  contact 
distance  d  is  not  exactly  the  sum  of  the  grain  radii  a  and  b.  If  the 

macroscopic  particles  have  a  magnetic  moment,  the  coagulation  process  Is 
considerably  enhanced  by  the  drift  motion  of  one  dipole  ~p^  in  the 

attractive.  Inhomogeneous  magnetic  field  Bg  of  the  other  dipole  ~pq. 
Experiments  Indicate  that  an  external,  homogeneous  magnetic  field  B0  directs 
the  agglomeration  of  particles  into  conglomerates  of  the  form  of  long  chains 
parallel  toUJ,.5*6 

In  order  to  provide  an  understanding  of  the  coagulation  of  ferromagnetic 
grains  in  dilute  suspension,  the  boundary-value  problem  for  the  spatial 

distribution  n(r)  of  similar  "field”  dipoles  ”p^  of  radius  r  a  in  the 

magnetic  field  Bg(r)  of  a  fixed  "sink"  dipole  ~pg  of  radius  r  •  b  shall  be 
analyzed,  when  a  homogeneous  magnetic  field  llj,  is  present.  The  latter  is 
assumed  to  be  parallel  to  the  z-axis  (Fig.  i)  and  is  in  spherical  coordinates 
(r,9,t)  given  by 

♦ 

B0  ■  B0  (cos0,  -  sinB,  0}  .  (45) 

The  sink  dipole  "pg  Is  taken  to  be  at  the  origin  T  ■  0  and  is  to  be  aligned 
with  T0  (Fig.  1),  so  that  in  spherical  coordinates 


The  magnetic  field  8a  »  -7(u  d  •r/4;rr 3)  of  the  sink  dipole  p„ 

D  O  8  j 

is  by  Eq.  (46) 

Bg  -  Ag  r“3{2cos0,  sine,  0 } ,  Ag  =  uoPg/4*  .  (47) 

The  Intensity  of  the  total  magnetic  field  B(r)  a  ~B0  +  "Bg(r)  at  the 

position  7  from  the  origin  (Fig.  1)  is  then 


B(r ,0)  -  (B2  +  2B  A  r" 3(3cos 20  -  1) 
o  OB 

(ABr"3)2(3cos23  +  D]  l/2,  d  <  r  <  -  .  (48) 

For  saturated,  ferromagnetic  grains  with  quasi-lnstantaneous  alignment12  in 
the  local  magnetic  field  B(r),  the  magnetic  dipole  force  is  derivable  from  the 

potential 17  ~Pa*B*> 

♦  ♦  ♦  ♦  ♦ 

pa  •  7B  -  7(pa.B),  pa  -  paB/B  .  (49) 


For  Isothermal  systems,  it  is  convenient  to  Introduce  the  dimensionless 
potential 


with 


*(r,0)  -  (paB(r , 0)  -  paB0]/KT 

(50) 

<»(r,e)  +  0,  r  ♦  •,  0  £  0  £  *  , 

(51) 

pa*7B/KT  -  7$  . 

(52) 

The  carrier  medium  in  which  the  magnetic  grains  are  suspended  is  assumed 
to  be  at  rest,  <v^>  ■  0.  According  to  Eqs.  (37),  (50),  and  (52),  the 
distribution  n(r,0)  of  field  dipoles  lTa  in  the  surrounding  d  <  r  <  •  of  the 
sink  dipole  7*}  (Fig.  1),  as  modified  by  thermal  diffusion  «7n,  magnetic 
field  drift  «7<i,  and  coagulation  at  the  contact  sphere  r*d*a  +  bls,  in 
the  stationary  state,  determined  by  the  elliptic  boundary-value  problem: 
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7*(7n  -  n7i|>)  *  0,  d  <  r  <  «, 


(53) 


n(r*d,9)*0,  O£0£ir,  (54) 

n(r  ■  »,0)  ■  n0,  0  £  0  <  ir,  (55) 

where  n  -  n(r,9)  due  to  azimuthal  ($)  symmetry,  and  Oq  Is  the  grain  density 
at  large  distances  r  »  d  from  the  sink  dipole  ‘p’ .  For  ij/(r,0)  =  0, 

D 

Eqs.  ( 53 )— ( 55 )  reduce  to  Smoluchowski ' s  boundary-value  problem,16  which  has 
the  simple  solution  n(r,6)  -  1^(1  -  d/r). 

The  magnetic  field  drift  term  «7iJ)  and  the  complex  r,  0-dependence  (48)  of 
the  potential  function  $(r,9)  render  the  boundary-value  problem  (53)-(55) 
nontrivial.  The  ansatz 

n(r ,9)  -  N(r,e)e*(r’0)  (56) 

reduces  Eq.  (53)  to  the  Laplace  equation  with  a  variable  coefficient  exp(i^), 

7»(e^7N)  -  0  .  (57) 

The  first  Integral  of  Eq.  (57)  is  proportional  to  the  grain  flux  X  * 
-D  exp(t(>)7N  -  -D(7n  -  n7i|i).  Hence, 

e*7N  -  {(3S/39)/r2sin9,  -( 3S/ 3r)/rsin9,  0}  .  (58) 


can  be  derived  from  a  stream  function  S  ■  S(r,9)  so  that  7*j  ■  0  and  Eq.  (57) 
is  satisfied.  Equation  (58)  yields 


(59) 


„  e  V  3 111  as  H  as 

72n  .  _  -s -  (—  —  -  —  — ) 

r zsin9  '‘ar  39  39  3r' 

The  right  side  of  Eq.  (59)  suggests  that  the  stream  function  Is  of  the  form 


S(r,9) 


fl(r,  9)e 


<Jt(  r ,  9 ) 


(60) 


which  eliminates  the  variable  coefficient  exp(-i|i)  from  Eq.  (59), 


V2N 


1 

r2sin9 


i£  12 

(3r  39 


ii  12 

39  3r5 


(61) 


In  order  to  homogenize  both  boundary  conditions  (54)  and  (55),  N(r,9)  Is 
decomposed  as 


N(r,9)  -  nQ(l  -  d/r)  +  g(r,9)  .  (62) 


By  Eqs.  (53)-(55),  (56),  (61),  and  (62),  the  function  g(r,9)  is  determined  by 
the  boundary-value  problem  for  a  Poisson  equation  with  homogeneous  boundary 
conditions: 


where 


72g  -  Q(r,9),  d  <  r  <  •*, 

g(r  -  d,9)  -  0,  0  1  ®  <.  *» 

g(r  -  •, 9)  •  0,  0  ^  9  <  w. 


Q(r,9) 


i  3j£  an 

r2sln9  ^  3r  36 


14  12 

ae  ar) 


(63) 

(64) 

(65) 


(66) 


In  the  field-free  case,  the  stream  function  is  S(r,9)  •  -  n^dcosa  for 
i|i(r,9)  ■  0.  Hence,  for  d»( r . 9 )  *  0, 


S(r ,9)  -  - 


<»(r  ,9) 

n  d  e  cos  9, 


n(r,9)  - 


-  n  dcos9, 
o 


(67) 


and 
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g(r, 9)  •  /  f  (Xr’.O'Wr.Ojr'.S'^irr'2  sine'dr'de'  .  (73) 

d  0 


Thus,  wc  obtain  from  Eqs.  (56),  (62),  and  (69)  the  following  analytical 
solution  for  the  spatial  distribution  n(r,6)  of  the  field  dipoles  ~pa  in  the 
superimposed  magnetic  field  Bg(r,9)  of  the  sink  dipole  'pg  and  the 
external,  homogeneous  field  BQ: 


n(r,9)  -  (n  (l-d/r)+  /  /  Q(r  ’ ,0  ' )G(r,9; r ' ,0 ' )2irr'  2sin9 ’dr ’d9 '  )e^  ’  ^  .  (74) 
0  d  0 
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This  result  shows  that  the  distribution  n(r,9)  of  ferromagnetic  grains  is 
composed  of  the  field-free  solution  multiplied  by  .  the  Boltzmann  factor, 
n0(l-d/r)exp(ii>),  and  a  source  solution,  g(r,  9)exp(  ij>),  which  vanishes  at  both 
boundaries  r  *  d  and  r  * 

The  current  density  j  ■  -D(7n  -  nV4>)  of  the  magnetic  grains  a  streaming 
to  the  interaction  sphere  r  *  d  around  the  sink  dipole  6  is  by  Eq.  (74), 


j  (r,  9)  -  -D(n  dr_2e  +7  ("(^(r ' ,  9  '  )G(r,  9 ;  r ' ,  9  ' )  2irr  ’  2sin9  'dr  fd9  •  le^  C  ’  9) .  (75) 
°  r  d  0 


As  an  illustration.  Figs.  2,  3,  and  4  exhibit  the  dimensionless  grain 
density  N(p,9)  ■  n(r,9)/n0  versus  the  dimensionless  distance  p  »*  r/d  1 

(from  the  source  dipole  8)  for  the  angles  0  »  0;  ir/4;  ir/2,  with  e0  = 
paB0/KT  *  1  and  e  5  paBg/KT  *  0.1;  l;  10  as  Interaction  parameters 

(&g  »  UoPB^1’*13^  Ic  18  seen  ttiat  tf,e  field  dipoles  a  are  distributed 

anisotropically  around  the  source  dipole  B»  caused  by  its  magnetic  self-field 
TTg(r,9)  and  the  external  (homogenous)  magnetic  field  B0,  where  X(o, 0)  > 
N(p**/4)  >  N(p,ir/2).  This  anisotrophy  effect  increases  strongly  with 

increasing  e  «  Bg  (Figs.  2,  3,  4),  but  Increases  only  slightly  with 

increasing  eQ  «  B0  in  the  interval  0. 1  <  e0  <  10  (for  the  latter  reason, 
only  e0  •  1  shown). 

Comparison  of  Fig.  4  with  Figs.  2  and  3  Indicates  that  the  magnetic  grain 
distribution  N(p, 9)  changes  qualitatively  for  large  interaction  parameters  e  > 
10.  The  distributions  N(p,9)  in  Figs.  2  and  3  are  of  the  "diffusion  type," 
whereas  the  distribution  N(p,9)  in  Fig.  4  is  controlled  by  the  (ac-6)  dipole- 
dipole  interaction.  Since  for  t  >  10  the  energy  paBg  of  the  field  dipoles 
a  in  the  magnetic  field  Bg  of  the  source  dipole  8  is  much  larger  than  the 
thermal  energy  KT,  a  large  amassment  of  field  dipoles  a  results  in  front  of 


the  reference  dipole  3  (Fig.  4).  In  view  of  the  binary  Interaction 
approximation,  Fig.  4  is  strictly  applicable  only  to  absolute  a-dipole 
densities  n(r,9)  for  which  the  average  binary  (a-a')  interaction  parameter 
eoa«  is  small  [Eq.  (7)]. 

In  view  of  these  observations,  it  is  to  be  expected  that  the  coagulation 
rate  yag  of  grains  with  magnetic  moments  p“a  and  ‘pg  depends  (i)  strongly 
on  the  magnetic  dipole-dipole  interaction  paramters  e  and  (ii)  on  the 
interaction  parameter  e0  of  the  dipoles  with  the  external  magnetic  field 
B0.  In  the  presence  of  a  strong  B0-field  the  coagulation  rate  Yag 
obviously  decrease  with  increasing  e  *  u0papg/4xd 3KT  so  that  grains  with 
large  magnetic  moments  pa  and  pg  coagulate  at  a  slower  rate  than  grains 
with  smaller  magnetic  moments  (in  agreement  with  experiments). 


COAGULATION  IN  MAGNETIC  FIELD 


The  local  density  nao  of  particles  (grains,  .chains)  with  magnetic 
moment  pa  decreases  due  to  coagulation  with  particles  of  magnetic  moment 


pg  and  local  density  ngQ  in  accordance  with  the  rate  equation 


16 


3n  /3t  =  -  T  y  „n  n  . 
ao  -  a6  ao  So 

0 


(76) 


The  summation  extends  over  all  particle  components  8  ^  a  including  3  =  a. 
Equation  (76)  defines  the  biliary  coagulation  coefficients  yag[ra3/sec] . 

For  the  analytical  calculation  of  the  coagulation  coefficients  yag  use 
is  made  of  the  numerical  result  that  the  spatial  distribution  na(r,9)  of 
a-dipoles  around  a  8-dipole  is  given,  in  good  approximation,  by  Boltzmann 
statistics  as  (dag  *  contact  sphere  radius  for  a-8  coagulation) 


n  (r, 0)  -  n  (1  -  d  Jr)  exp[<J»  (r,0)] 
a  ao  an  an 


(77) 


where 


Wr’9) 


p  B  p  B 
a  o  .  a  o 


KT  KT 


+  -  —  )[1  +  3(2  -1  +  l)(-2.  -  I)"2  cos20]l/2  (78) 


and 


rs  5  (*S/B0>l/S  ■  (“opS/4*,o,l/3  •  (7,) 

Equation  (77)  follows  from  the  complete  solution  (74)  by  neglectlon  of  the 
g(r,0)  contribution  (78)  with  vanishing  boundary  conditions  (Eqs.  ( 63 ) — ( 66 ) ] . 

-IT  n 

In  this  approximation,  the  flux  tag  -  2nf  jr t  agr^sinSdS  of  a-dipoles 


to  the  8-dlpole  is 


*«8(r)  *  '  Da2imaoda8  fW  eX?  f^ft(r*9)]  sin0  d9 


since  Che  a-dipole  current  density  is  jag  *  -  Da(7na  -  » 

-  DanaodaSr~  exp(  (iiag)er  by  Etl*  (77).  The  number  of  a-dipoles  which 


reach,  per  unit  time,  the  contact  sphere  r  *  dQg  of  one  8~dipole  is 


°oo  •'as  ‘  •  Ca3VC  •  da0>  ’ 


L,  a  «  8;  C  "*  2,  a  *  8.  (81) 

aB 


Cag  is  an  accommodation  coefficient,  which  considers  that  not  only  the 


a-dipole  but  also  the  8-dipole  is  in  thermal  motion.  16  Equations  (80)  and 


(81)  yield  for  the  coagulation  coefficient  in  an  external  (homogeneous) 


magnetic  field  the  formula: 


p  B  1  p  B  (2p3  +1)  . 

'’.a  *  /  "Pl-prO  -  »;J)f'+3  :  »2]l/2td»  (82) 


where 


°.8  3  da8/c8  ’  <4’daV,o/“oPB) 


and  a  -  cosB.  For  the  practical  evaluation  of  the  integral  in  Eq.  (82)  it  is 


noted  that  the  dimensionless  parameter  pag3  is  very  large  even  for  moderate 
external  B0-fields,  e.g.,  pag3  —  104  for  dag  —  10_7m,  pg  —  1018Am2, 

and  B0  ■  1  Vsec/m2.  Accordingly, 


p  B  n 


Ya0  "  4*Ca8d<*0Da  exp("  n  j  ^  *3...(2n+l)  (3  *  °a8  >  1 


lipp,  »  UPP- 

~  _  .  n  r  o  a  8  >r.  r  (—6) _ f  o*a 

YaS  *  «B  a8  a  4„d3  KT)[  Jjl  0...(2n+l)  ^  3  >  ^ 

aB  aB 


V  Vs /4’daS 


WWW 


Although  Eq.  (82)  predicts  a  B0-dependence  in  general,  the  coagulation 

coefficient  yag  is  no  longer  B0-dependent  in  the  limit  of  large 

B0-fields,  Pag3  >  1  or  B0  >  UoPe/^^aB 3*  Equation  (85)  shows  that 

Yag  decreases  exponentially  with  increasing  ratio  of  dipole  interaction 

energy  UoPaPS^^aS 3  to  thermal  energy  KT.  Hence,  grains  or  chains 

with  large  dipole  moments  coagulate  slower  than  those  with  small  dipole 
moments,  in  the  presence  of  a  strong  external  magnetic  field  B0  (Eq.  (85)1. 

In  the  absence  of  an  external  magnetic  field,  T}0  •  0,  we  have  paB  =  0 
by  Eq.  (83).  In  this  case,  Eq.  (82)  gives  for  the  coagulation  coefficient 

1  exp +  3o2)l/2]do  ,  B  =>  0  .  (86) 

4ird3  KT  ° 

oB 


Ya8  -  4,fCaBda8Da  [ 


Since  0  <  3o2  <  3  for  0  <_  o  <_  1,  the  integral  in  Eq.  (86)  can  no  longer  be 
■i  approximated  by  a  rapidly  converging  series.  For  particles  without  dipole 

interactions  (rigid  spheres),  0o>g  *  0,  Eq.  (86)  reduces  to  Smoluchowski ' s 
formula  ya8  *  4irCagdagDa. 1 6 

The  integral  functional  in  Eq.  (86)  indicates  that  the  coagulation  rate 
ya8  increases  essentially  exponentially  with  Increasing  ratio  of  dipole 
interaction  energy  u0papg/4ird0g3  to  thermal  energy  KT,  if  an  external 
magnetic  field  is  not  present.  Thus,  in  the  case  B0  -  0,  ya8  is  larger 
than  in  the  presence  of  a  strong  external  magnetic  field,  B0  » 
WoP8^ff<*aB3,  *n  t*ne  latter  case,  y0g  decreases  exponentially  with  the 
same  energy  ratio  (Eq.  (85)]. 

The  stabilizing  effect  of  an  external  magnetic  field  B0  concerning  the 
decay  of  a  magnetic  aerosol  by  coagulation  is  experimentally  established.  In 
particular,  homogeneous  magnetic  fields  B0  _>  10~2  Vsec/m  are  applied  to  the 
chambers  of  ovens  for  the  generation  of  ferromagnetic  aerosols,  in  order  to 


reduce  Che  coagulacion  rate  and  to  align  the  ferromagnetic  colloid  chains.  j*y 
Eq.  (76)  the  relaxation  time  Tac  for  coagulatlve  decay  of  the'  density 
nao  of  dipole  particles  with  magnetic  moment  pa  is 


c 

T 

a 


-1 


aS 


UV«  0 


(87) 


Another  method  of  reducing  the  decay  rate  of  ferromagnetic  aerosols  is  to 
spray  the  ferromagnetic  particles  with  electric  charges.  The  resulting 
Couloumb  repulsion  of  the  magnetic  particles  effectively  reduces  coagulation 
as  has  been  demonstrated  experimentally. 17 

The  presented  theory  should  be  considered  as  a  first  step  towards  a 
quantitative  understanding  of  the  coagulation  of  magnetic  aerosols  in  external 
magnetic  fields.'  At  higher  aerosol  densities,  many  -  dipole  interactions  have 
to  be  considered  (in  addition  to  the  binary  dipole  interactions)  in  the 
analysis  of  the  coagulation  coefficients.  In  principle,  this  can  be 
accomplished  by  means  of  the  kinetic  equations  ( 2 1 ) — ( 24 )  or  the  corresponding 
macroscopic  transport  equations  (37 )— ( 42 ) ,  with  self-consistent  magnetic 
dipole  interactions. ' 
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FIG.  2:  Dimensionless  grain  density  N(p* 9)  versus  p  ■  r/d  l  for  9 
ir/2,  e„  »  1,  and  e  ■  0. 1. 


FIG.  3:  Dimensionless  grain  density  N(p»9)  versus  p  -  r/d  £  I  for  9 
t/2,  e0  -  1,  and  e  ■  l. 


FIGURE  CAPTIONS 

FIG.  1:  Field  dipole  pa  in  the  magnetic  field  B(r,0)  of  a  source  dipole 

pg  and  a  homogeneous  magnetic  field  B0. 

FIG.  2:  Dimensionless  grain  density  NCp,  9)  versus  p  =■  r/d  >  1  for  0  **  0;  n/ 4; 
it/2,  eQ  *  1,  and  e  =»  0.1. 

FIG.  3:  Dimensionless  grain  density  N(o,  0)  versus  p  *  r/d  1  for  0*0;  ir/4; 
n/2,  e0  -  1,  and  e  -  1. 


FIG.  4:  Dimensionless  grain  density X(p. 9)  versus  p  -  r/d  ^  1  for  9*0;  tt/4; 
ir/2,  eQ  *  1,  and  e  *  10. 
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STATISTICAL  DISTRIBUTION  OF  RANDOM  ELECTRIC  FIELDS  IN  CHARGED  PARTICLE  GASES 

H.E.  Wilhelm 
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ABSTRACT 

The  temperature  (T)  and  density  (n)  dependent  probability 
distribution  W  *  W(E;T,n)  of  the  collective  electric  fields  E  in  an 
ideal  gas  of  charged  particles  each  carrying  the  same  charge  £  (electrons: 
£■-^<0;  ions:  £ ■  Z£q  ^  0)  is  calculated  from  first  principles  of  stati¬ 
stical  mechanics.  It  is  shown  that  (i)  the  collective  electric  fields 

act  at  distances  larger  than  the  characteristic  repulsion  distance  D  * 

2  1/2 

(KT/4Trn£  )  of  like  charged  particles,  and  (ii)  the  average  collective 

1/2 

electric  field  is  *  (12irnKT)  for  ideal  gases  of  particles  of  the 
same  charge  £.  Thus,  in  a  thermal  equilibrium  gas  of  like  charged  par¬ 
ticles,  the  longitudinal  microfields  are  considerably  stronger  than  in  a 
quasi-neutral  plasma,  since  in  the  latter  the  random  field  effects  of  the 
negative  electrons  and  positive  ions  nearly  compensate  each  other.  Finally, 
the  interrelation  between  (average)  kinetic,  interaction,  collective  field, 
and  electric  self  energies  is  discussed  for  charged  particle  gases. 


'men. v 


INTRODUCTION 


The  probability  distribution  of  the  stochastic  electric  fields  E 


produced  by  charges  e  5  0  of  one  kind  in  random  thermal  motion  is  of 


considerable  interest  since  one -component  charged  particle  gases  such  as 


electron  or  ion  gases  are  employed  in  many  technical  applications  of  physi¬ 


cal  electronics  (von  Ollendorff  1957).  In  a  macroscopically  homogeneous 


charged  particle  gas  of  density  n  and  temperature  T,  collective  particle 


interactions  occur  at  distances  r  >  D,  since  the  minimum  wave  length  of  the 


random,  thermally  excited  longitudinal  charge  waves  is  of  the  order  X  -  D 

m 


of  the  characteristic  repulsion  distance  for  particles  with  the  same  charge  e. 


D  -  (KT/4irne2)1/2 


By  means  of  Poisson's  equation,  the  random  collective  field  amplitudes  in 


an  arbitrary  direction  "i"  can  be  estimated  from  the  random  particle  density  n  as 


E^/D  *  ±4iren 


Equation  (1)  and  (2)  show  that  an  equipartie-Lon  between  random  collective 


field  and  thermal  energies  exists  on  the  average  in  a  gas  of  like  charges. 


<E*/8*>  -  nKT/2 


A  rigorous  derivation  of  Eq.  (3)  based  on  the  Markov  method  has  been  given  by 


Mints  (1957)  for  a  gas  of  electrons  in  thermal  equilibrium. 


By  the  fundamental  axiom  of  statistical  mechanics  of  ideal  systems  in 


thermal  equilibrium,  all  equilibrium  distributions  can  be  derived  without 


consideration  of  the  interactions  which  bring  about  the  equilibrium  (Tolman 


1938).  By  extending  this  principle  for  many-particle  systems  with  discrete 


energies  to  continuous  media  with  random  energy  densities  u  »  E(r,t)"’/8Tt, 


we  derive  the  probability  distribution  W(E)  of  the  collective  electric  fields 
E(r,t)  in  charged  particle  gases.  These  are  assumed  to  be  "ideal",  i.e.,  the 
Coulomb  repulsion  energy  &  / r  is  assumed  to  be  small  compared  with  the 
thermal  energy  KT, 

Y  «  e2n1/3/KT  =  1.670*10" 3Z2n1/3T_1  <<  1  .  (4) 

Among  the  results  reported,  it  is  shown  that  the  longitudinal  microfields 
in  gases  of  like  charged  particles  are  much  . larger  than  those  of  quasi-neutral 
plasmas,  in  which  the  random  electric  fields  of  the  negative  electrons  and 
positive  ions  nearly  compensate  each  other. 


PROBABILITY  DISTRIBUTION 


Subject  of  the  considerations  is  a  homogeneous  gas  of  volume  0  containing  n 
charges  z  per  unit  volume.  In  thermal  equilibrium,  the  kinetic  energy  density 
of  the  charges  of  mass  m  and  velocities  v^ is  given  by 

N  ,  - 

<  l  j  mv  >  -  |  nP.KT  ,  N  *  n.Q  (5) 

y-1 

During  the  random  thermal  motions  of  the  charged  particles,  a  continuous  trans¬ 
formation  of  kinetic  particle  energy  into  potential  electric  energy  occurs,  and 

vice  versa,  due  to  the  particle  interactiond  through  their  longitudinal  Coulomb 

2 

fields  (transverse  electromagnetic  interactions  are  negligible  for  me  >>  KT) . 

As  has  been  shown  -first  by  Mints  (1957),, equipartition  of  average  random  electric 
and  kinetic  energies  exists  in  statistical  equilibrium  [for  a  thermodynamic  deri¬ 
vation,  see  Eq.  (46)]  of  a  gas  of  particles  with  the  same  charge, 

<t2/8ir>  -  |  nKT  .  (6) 

The  electric  field  E(r,t)  produced  collectively  by  the  charges  at  any  point 

4  >f  2 

r  e  D  and  the  field  energy  density  u  »  E(r,t)  /8ir  fluctuate  with  time  t  about  the 

average  values  <fe>  -  t)  and  <d2/8it>  i*  0,  respectively.  The  proposed  problem  is  to 

derive  the  probability  wdjd3!  for  finding  the  collective  field  fluctuation  E  in 

the  volume  element  d3!  -  dE  dE  dE  about  the  point  1  ■  (E  ,E  ,E  )  of  the  field 

x  y  z  x  y  z 

space  subject  to  the  thermal  equilibrium  conditions  (5)  and  (6). 

In  order  to  determine  experimentally  the  collective  microfield  distribution 
W(l)  -  W(l2/8ir)  in  a  homogeneous  and  isotropic  gas  of  charged  particles,  one  would 
have  to  measure  the  fluctuating  field  E  or  the  fluctuating  energy  density  E  / 8tt  in 
the  vicinity  6  r  of  a  fixed  field  point  reft  at  consecutive  times  ty  *  v@^. 


5 


v  »  1,2,3, .. .V  within  experimental  errors  At  <<  G  ,  where  0  is  a  time  interval 

V  V  * 

which  is  large  compared  with  the  correlation  time  of  the  stochastic  field  so  that 

<l(r,t)l(r,t  +  0y)>  “0  (within  these  limitations,  the  magnitude  of  0^  may  be 

changed  from  one  measurement  v  to  the  next  v  +  1) .  In  a  large  number  of  such 

-*2 

measurements,  N  ■*  ®,  the  energy  density  E^/8ir  would  be  observed  times,  ...» 

2  2 
the  energy  density  E  /Sir  would  be  observed  N  times,  etc.,  where  1  /8n  means  an 

a  a  a 

experimental  value  measured  with  an  error  A(£2/8it).  The  resulting  step-shaped 

energy  distribution  N  »  N  (E2/8it)  is  represented  by  the  partition 

a  qi  ci 


H, 


*•« 


N 

a 


N, 


where 


tj/8ir  l2 /8ir  l2/8n  ...  l2/8ir  . . .  t£/ 


N1  +  N2  +  N3  +  *  *  *  Na+  •  •  •  Njj  *  R 


M 

8x 


(7) 


(8) 


Hjtj/ftr  +  N2^/8ir  +  N3I3/8X  +  . . .  Nl2/8ir  +  ...  N^/8tr  -  IV<E2/8ir>.  (9) 

N  is  the  total  number  of  measurements  ( N  •*  «•)  and  #<E2/8tt>  is  the  total  field 
energy  density  measured  in  the  N  independent  observations.  The  entire  energy 
density  N<%2 / 8ir>  can  be  distributed  in  a  large  number  II  of  ways  over  sets  {N  } 


aN 


of  numbers  If  .  By  elementary  combinatorics  (Tolman  1938) , 


H  -  N\  /  Nj^I  N2!  N3!  ...  N#!  ...  Njj!  .  (10) 

The  energy  distribution  N  (!2/8w)  observed  in  statistical  equilibrium  is  the 


a  a 
tZ 


most 


probable  one.  Thus,  Na(l^/8ir)  is  determined  by  the  condition  for  a  maximum  of  i) 
the  number  H  of  realizations  or  ii)  the  entropy  S  -  lnH,  subject  to  the  constraints 
(  8 )  and  (  9 ) . 


6 


Accordingly,  we  determine  the  probability  distribution  from  the 


maximum  of  the  function  lnll  =  f(N^), 


M 


Inn  -  N(lnN  -  1)  -  Y  N  (InN  -  1) 

L.  a  a 
a*l 


with 


M 


I ».  - »  • 


a*l 

M 


N  -*■  « 


y  N  1  /8tt  -  iV4  nKT 
a  a  2  . 

a*l 


0 


(ID 

,  (12) 

,  (13) 


as  constraints.  Eq.  (12)  holds  by  definition  of  N,  whereas  Eq.  (13)  holds  for  a 
large  number  N  of  measurements  and  the  average  energy  density  <E  /8ir>  of  Eq.  (6). 
Addition  of  the  constraints  (12)  and  (13)  multiplied  by  the  Lagrangian  multipliers 
and  -y  to  Eq.  (ll)  leads  to  the  compact  maximum  conditions  for  Inn, 

3F(N  )/3N  -  0  ,  32F(N  )/3N2  <  0  ,  o  -  1.2....M  ,  (14) 

a  a  a  a 

where 

M  MM 

F(N  )  -  HlnN  -  1)  -  7  N  (InN  -  1)  -  A  Y  N  -y  F  N  l2/8ir  .  (15) 

(X  -  Qt  (X  y  (X  “  _  (X  (X 

a*l  a*l  a*l 


The  solution  of  Eq.  (14)  gives  the  distribution  N  of  the  "discrete"  energy  densi- 

a 

ties  l2/8ir  in  the  form 
a 

-yl2/8ir 


N 


Ae 


A  -  e-(1+A> 


(16) 


Henceforth,  the  subscript  a  is  dropped  since  can  be  any  point  1  in  the  field 
space.  The  dimensional  constants  A(A)  and  y  are  then  given  by  the  normalization 
conditions  (12)  and  (13) , 


.  r°°  -yt  /8ir.  _2  „ 

A  J  e  4wE  dE  ■  N 

0 


1 

I 

I 


(17) 


7 


/  (E2/87r)e“y^  /8lT4irE2dE  =  <V~ 


,  CIS) 


A  »  (8TT2nKT)‘3/2iV 


V  =  1/nKT 


For  this  normalization,  which  still  contains  the  number  N  of  measurements,  the 
probability  distribution  (16)  for  the  raicrofield  energy  density  is 


W  (£2/8ir)  *  T  — — rrr-  e 

N  (sAkt)372 


-£2/  8TmKT 


•  (20) 


In  theoretical  applications,  one  is  interested  in  the  probability  dP  =*  W(E)d  E 
for  finding  a  microfield  £  in  the  volume  element  d3£  about  the  point  E  of  the  field 
space,  with  the  normalization  /dP  -  1.  The  corresponding  distribution  function 
W(£)  of  the  collective  microfield  £  is  obtained  by  renormalization  (N  -*■  1): 


W(£)  -  (8n2nKT)'3/2e~^  /8mlKT 


.  (21) 


This  fundamental  distribution  has  the  form  of  a  Gaussian,  l.e.,  all  its  moments 
exist,  e.g.. 


<£°>  -  ///£°W(£)  d3£  -  £ 

<£1>  -///* £1w(£)d3£  -  £ 

—  00 

<£2>  -  ///£2w(£)d3£  -  12  tokt 


,  (22) 


,  (23) 


•  (24) 


The  most  probable  (E^)  and  the  r.m.s.  (Ew)  collective  microfields  are  by  Eqs.  (21) 


and  (24) 


Ejj  -  (8imKT)- 


,  (25) 


Ey  -  (12imKT) 


(26) 


Far  considerations  concerning  the  fluctuation  of  the  collective  microfifeld 


E(t)  at  a  point  reft  with  time  t,  temporal  averages  can  be  defined  by 

— — -  .  +T 

|E|  -  lia  £/  |l(t)|dt  ,  (27) 

T-*“  -T 

t1  *  lim  j-  /l(t)2dt  .  (28) 

T-voo  _-f 

The  fluctuation  of  £(t)  is  defined  by  AE(t)  =  E(t)  -  l(t)  with  E(t)  =  5.  In 
stationary  equilibrium,  the  time  averages  are  identical  with  the  ensemble  averages 
By  Eq.  (21),  the  mean  square  (temporal)  fluctuation  of  £(t)  is 

AE2  -  I2  -  |£|  -  (3  -  ®)4irnKT  .  (29) 

2  1/2 

TABLE  I  compares  the  r.m.s.  field  and  the  r.m.s.  fluctuation  (AE  )  of 

2/ 

the  collective  microfield  with  the  nearest  neighbor  Coulomb  field  Eq  *  2ir(4/5) 
2/3 

|e|n  in  dependence  of  the  particle  density  n  for  an  ideal  electron  gas 
4  2  1/2 

(y  <<  1,  T  ■  10  °K).  It  is  seen  that  E^  and  (AE  )  are  one  to  two  orders  of 

18  -3 

magnitude  larger  than  E^  in  the  range  of  ideal  gas  densities  n  <  10  cm  .  For 
these  reasons,  the  Coulomb  field  Eq  represents  a  small  contribution  to  the 
microfield  in  ideal  gases  of  charged  particles.  The  result  E^  >>  Eq  is  readily 
understandable  since  for  ideal  conditions 


E2/eJ  -  (ir/3)(4/15)4/3- 


e2n1/3 

KT 


<  y  <<  1 


(30) 


The  probability  for  observing  a  collective  microfield  with  intensity  E  *  | E | 

2  *► 

in  the  range  between  E  and  E  +  dE  is  P(E)dE  ■  W(E)4itE  dE,  where  W(E)  is  given  by 

-1  2  -1/ 

Eq.  (21).  The  maximum  of  the  probability  density  P(E)  is  P(E^)*4exp  (8n  nKT) 
by  Eq.  (25).  Accordingly,  the  normalized  probability  density  is  P(E)/P(EM)  * 
x  (E2 /  8imKT)exp[l  -  E2/87mKT]  <  1.  Figure  1  presents  P(E)/P(EM)  -  P(E2/8TinKT) 
versus  E  /8nntCT.  This  distribution  is  a  displaced  Gaussian  with  a  maximum 
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P(E^/8"nKT)  =  1  at  E  =  E^.  The  most  probable  electric  field  =  (8-rnKT) 

1/2 

increases  proportional  (nT)  with  increasing  nT  values.  The  increasing 
quantitative  importance  of  the  collective  microfield  E  in  ideal  charged  par¬ 
ticle  gases  with  higher  pressures  p  =  nKT  is  obvious. 

On  the  other  hand,  in  a  quasi-neutral  electron  (e)  -  ion  (i)  plasma 

3  -1 

(ng  =  Zn^  ,  the  microfield  energy  Up  differs  by  a  factor  of  order  (neD+) 

2 

from  the  microfield  energy  E^/.frr  [Eq .  (26)  ]  of  the  charged  particle  gas,  where 
2  2  1/2 

D  =  [KT/4ir(n  e  +  n.£. )1  4  D  is  the  Debye,  shielding  radius  (Debye  and 

±  e  e  l  l 

3  -1 

Hueckel  1923).  Since  (n  D  )  <<  i.  for  ideal  conditions,  the  microfields  in 

e  - 

plasmas  are  small  in  comparison  with  those  in  ideal  gases  of  like  charged 
particles.  In  a  plasma,  the  random  electric  fields  of  the  negative  electrons 
and  positive  ions  compensate  each  ether  nearly  completely.  This  incomplete 
statistical  compensation  is  the  physical  reason  why  0  <  Up  <<  E^/8»  . 


ENERGY  RELATIONS 


A  gas  of  charged  particles  in  thermal  equilibrium  at  a  temperature  T  ex- 

3 

hibits  various  macroscopic  energies,  the  average  kinetic  energy  <K>  -  -^nlCTft,  the 
average  electric  field  energy  <U>  ■  <E^/8tt>0,  the  average  interaction  energy  <4» 
and  self energy  <*K>  of  the  charged  particles.  In  order  to  derive  the  interrela¬ 
tion  between  these  energies,  the  formation  of  the  gas  by  an  electric  charging 
process  is  considered.  For  this  purpose,  we  assume  that  the  charged  particles 
are  initially  dispersed  at  infinity  where  they  have  only  self energies.  The  gas 
is  then  built  up  by  moving  one  charge  after  the  other  from  infinity  into  the 
volume  0,  which  requires  work  against  the  resulting  Coulomb  field  of  the  charges 
already  present  in  ft.  The  thermodynamics  of  the  charging  process  is  illustrated 
by  i)  a  reversible  isothermal  and  ii)  an  adiabatic  or  isotropic  model. 

The  electric  charging  work  expanded  in  moving  N  charges  e  against  their 
collective  Coulomb  field  from  infinity  into  the  (finite)  volume  ft  is  (*  designates 
exclusion  of  terms  with  p  ■  v) 


1  w 

il 


N  .  N 


p-1  v 


I  V 
1 


r  -r  , 

p  v ' 


i-l 


p,v 


e  ft 


(31) 


where  r^  (r^)  are  the  position  vectors  of  the  p-th  (v-th)  charge  <L  in  the  volume 
ft,  respectively.  The  collective  microfield  of  the  N  charges  at  a  field  point 
(r,t)  is  the  superposition 


l(r,t)  -  l  tCr.t)  (32) 

P-1 

where  E^(r,t)  are  the  individual  Coulomb  fields  produced  at  the  field  point 
(r,t)  by  the  p-th  charge.  By  Eq.  (32),  the  electric  field  energy  of  the  gas  ft  is 
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(33) 

(34) 

(35) 


are  the  (e-e)  interaction  energy  and  the  (e) self energy  of  the  charged  particle 
gas,  respectively.  Comparison  of  Eq.  (31)  with  Eq.  (33)  reveals  the  interre¬ 
lation 


U  -  ¥  -  4  *  A 


(3b) 


Thus,  we  see  that  the  field  energy  U  is  the  sum  of  the  interaction  energy  4  and 
the  self energy  [Eq.  (35)].  The  charging  work  A  leads  to  an  increase  of  the  in¬ 
teraction  part  4  of  the  field  energy  U  [Eq.  (36)].  The  selfenergy  ¥  of  the  charges 
is  Independent  of  the  spatial  locations  of  the  charges,  i.e.,  ¥  is  the  same  before 
and  after  the  charging  process. 

Another  independent  energy  relation  is  obtained  by  multiplication  of  the 

2+  2 

coupled  Newtonian  equations  for  the  accelerations  d  r  (t)/dt  of  the  u-th  charge 

V  •  v 

and  the  v-th  charge  by  their  respective  velocities  v^  y(t)  •  dr^  y(t)/dt  and  sub- 

I 

sequent  summation  over  all  particles  y  and  v.  The  resulting  expression  can  be 
brought  into  the  form  d(K  +  4)/dt  *  0,  which  demonstrates  that  the  sum  of  kinetic 
( K )  and  interaction  (4)  energies  is  an  invariant  Ho, 

K  +  4  -  Hq  (37) 

where 

N  l 

K-  I  f»I2  (38) 

U-l 

and  4  ■  A  Is  defined  in  Eq.  (31).  Eq.  (37)  expresses  the  conservation  of  kinetic 
K  and  interaction  4  energies  in  a  gas  of  charged  particles  which  interact  by  longi¬ 


tudinal  Coulomb  fields 
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The  thermodynamic  functions  of  the  gas  depend  in  general  on  the  volume  3, 
the  number  N  of  particles  in  0,  and  the  particle  averages  of  the  random  kinetic 
energies  and  the  random  field  energy  densities  E2/8tt.  Accordingly,  we 

assume  *  Ut^(T,e,N)  for  the  thermal  energy  and  S  *  S(T,e,N)  for  the  entropy, 
where 

3KT/2  *  <-| mv2>  ,  e  -  <E2/8tt>  .  (39) 

For  gas  formation  by  isothermal  reversible  charging,  the  volume  0  is 
embedded  into  a  heat  bath  of  temperature  T.  The  transfer  of  dN  charges  e  from 
infinity  into  the  cavity  n  requires  on  the  average  the  charging  work  dA  *  d<U  -  ¥> 
»  d<U>  by  Eq.  (36),  and  their  thermalization  at  a  temperature  T  consumes  on  the 
average  the  energy' dU  -  jKTdN.  The  difference  of  these  energies,  dQ,  is 
supplied  by  the  heat  bath.  In  accordance  with  the  first  law  of  thermodynamics 

dQ  -  dUCh  -  d<U>  (AO) 

since  no  other  than  electric  charging  work  is  performed  on  the  system  (dfi  =  0) . 

The  associated  entropy  dS  *  dQ/T  is  a  complete  differential, 

ds  ’  f  ^r(uth  *  <u>>dT  +  t  ^(uth  "  <l>)de  +  ^(uth  "  <U>)dN  •  <*i> 

Application  of  the  condition  3£3,j,S  *  3 ^,3 ^ S  to  Eq.  (41)  yields  the  partial  differ¬ 
ential  for  constant  N  and  T, 

3Uth/3e  -  3<U>/3e  .  (42) 

Since  =  0  for  e  »  0  (no  thermal  energy  in  0  before  charging) ,  the  integral 
of  Eq.  (42)  is 

th  .  ..  . 
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A 
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Eq.  (43)  could  have  been  derived  by  other  thermodynamic  gas  formation 
processes,  e.g.,  by  adiabatic  charging  of  the  cavity  0.  In  this  case  dQ  =  0, 
and  by  Eq.  (40) 


dQ  *  dU  -  d<U>  -  0  : 


<U>  =*  U 


Finally,  <U>  can  also  be  determined  as  that  equilibrium  value  which  maximizes 


the  entropy. 


dS  -  T'1[dUth  -  d<U>  ]  *  0  : 


<U>  *  U 


.  (45) 


Eqs.  (43)  -  (45)  indicate  that  an  equipartition  between  thermal  energy  and 
average  microfield  energy  exists  in  statistical  equilibrium.  This  fundamental 
result  is  explicitly  (Mints  1957) . 

|  NKT  -  0<E2/8it>  .  ( l 


CONCLUSION 


In  ideal  gases  of  charged  parcicles,  Che  distribution  function  of  the 
collective  microfields  is  strongly  temperature  and  density  dependent.  For 
typical  temperatures  and  densities  of  ideal  systems,  the  r.m.s.  collective 
microfield  is  by  orders  of  magnitude  larger  than  the  quasi-static  Coulomb 
field.  In  statistical  equilibrium,  a  balance  among  (average)  kinetic  particle 
and  collective  microfield  energies  exists,  which  is  independent  of  the  process 
of  the  formation  of  the  charged  particle  gas. 

The  derived  formulas  for  the  average  microfield  and  energy  are  appli¬ 
cable  to  ideal  gases  consisting  of  one  species  of  charged  particles.  Examples 
are  electron  gases  in  highpower  tubes  and  diodes,  and  non-neutralized  beams  and 
clouds  of  charged  particles  in  outer  space.  Our  results  are  not  applicable  *- 
plasmas,  since  in  these  the  electron  and  ion  components  have  no  independent 
existence  (coupling  of  the  negative  and  positive  charges  through  the  self- 
consistent  field  and  quasi-neutral  or  quasi-compensated  behavior) . 
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Abstract 

Based  on  the  Galilean  relativity  principle  and  Maxwell's  equations, 
electromagnetic  field  equations  are  derived  for  inertial  frames,  in  which  the 
substratum  of  the  electromagnetic  waves  flows  with  arbitrary  velocity  |w|  <  c 
(velocity  of  light) .  It  is  demonstrated  that  the  electromagnetic  field  equa¬ 
tions  with  electromagnetic  substratum  flow  are  strictly  covariant  against 
Galilei  transformations.  Wave  equations,  conservation  and  invariance  theorems, 
and  boundary  conditions  are  derived  for  the  electrodynamic  fields  in  presence 
of  electromagnetic  substratum  flow.  Initial-boundary-value  problems  are  solved 
for  electromagnetic  signal  propagation  and  induction  in  the  substratum  by  an 
Integral  equation  method.  Physical  effects -for  the  measurement  of  the  veloc¬ 
ity  field  of  the  electromagnetic  substratum  are  discussed.  Maxwell's  conception 
that  his  equations  refer  to  a  frame  of  reference  with  resting  electromagnetic 
substratum  is  confirmed,  and  it  is  shown  that  Maxwell's  equations  are  also 
applicable  to  inertial  frames  with  small  substratum  velocities,  |w|  <<  c. 
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INTRODUCTION 


Since  the  discovery  of  the  2.7°K  cosmic  microwave  radiation  [Wilson,  1980] 


Maxwell's  original  ideas  on  the  propagation  of  electromagnetic  waves  in  the  so-called 


electromagnetic  ether  have  become  of  renewed  interest.  According  to  Maxwell, 


Heaviside,  Lorentz,  and  Poincard  the  electromagnetic  field  equations  refer  to  a 


system  of  reference  in  which  the  carrier  of  the  electromagnetic  waves  is  at  rest 


(Whittaker,  1951].  By  comparison  with  ocher  wave  phenomena,  this  restriction  appears 


to  be  physically  and  mathematically  necessary,  since  the  Maxwell  equations  do  not 


contain  explicitly  the  velocity  field  w  of  the  electromagnetic  ether  [Lorentz, 


1909].  The  most  familiar  "material"  properties  of  the  ether  are  the  electric  (ec 


10“9/36x  Asec/Vm)  and  magnetic  (u0  »  4irxl0“7  Vsec/Am)  permeabilities,  the  wave 
resistance  ZQ  «  (uo/Co^1^2  *  376.731  0,  and  the  velocity  of  light  cQ  - 

(Uoeorl/2  “  3xl08m/sec  [Stratton,  1941].  The  difficulty  of  observing  the 


substratum  by  other  chan  electromagnetic  experiments  (e.g.,  measurements  of  the 


velocity  of  light  and  frequency  shifts)  is  probably  due  to  an  extremely  small  inter¬ 


action  cross  section  of  particles  with  the  substratum.  Evidence  for  this  is  given  by 


Flzeau's  [1851]  ether  drag  experiment  which  shows  that  the  ether  is  not  noticeably 


carried  along  by  liquid. matter  flowing  in  tubes. 


As  noted  by  cosmologlsts,  e.g.,  Mansouri  and  Sexl  [1977],  "The  validity  of  the 


principle  of  relativity  (which  assumes  that  the  velocity  of  one  and  the  same  light 


signal  has  the  same  value  vs  -  c0  in  all  (•)  inertial  frames)  seems  to  be  less 


evident  now  than,  say  twenty  years  ago.  The  discovery  of  the  cosmic  background  radi¬ 


ation  has  shown  that  cosmologlcally  a  preferred  system  of  reference  does  exist."  In 


critical  analyses,  Ives  [1938,  1948],  Builder  [1958a,  b],  and  Janossy  [1953,  1963] 


demonstrate  that  the  ether  effects  cancel  out  in  the  measurements  of  the  velocity  of 


light  by  Michelson-Morley  [1887],  Morley-Miller  [1905],  and  their  modern  versions 


[Jaseja  et  al.,  1964;  Vessot  et  al.,  1979]  (using  signals  sent  out  and  back).  On  the 


other  hand,  they  show  that  the  experiments  of  Sagnac  [1937]  and  Dufour-Prunlen  [1937] 
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(rotating  interferometers),  Ives-Stillwell  [1938],  and  Michelson-Gale  [1925]  support 
the  electromagnetic  ether  concept.  The  isotropic  microwave  background  radiation  in 
the  universe  appears  to  indicate  thermal  excitations  of  the  ether  at  a* nearly 
isotropic  temperature  Tg  -  2.7°K  [Wilson,  1980]  and  ether  velocities  in  the 
terrestrial  space  of  the  order  w  -  105m/sec  [Henry,  1971], 

The  impossibility  to  carry  the  denial  of  the  electromagnetic  substratum  in  the 
special  theory  of  relativity  over  into  the  general  theory  of  relativity  was  clearly 
recognized  by  Einstein  [1921]:  “According  to  the  general  theory  of  relativity  space 
without  ether  is  unthinkable;  for  in  such  space  there  not  only  would  be  no  propaga¬ 
tion  of  light,  but  also  no  possibility  of  existence  for  standards  of  space  and  time 
(measuring  rods  and  clocks),  nor  therefore  any  space-time  in  a  physical  sense." 
However,  Builder  [  1958a, b]  demonstrated  for  the  first  time,  by  means  of  quantitative 
theoretical  arguments,  that  relativity  theory  without  electromagnetic  ether  leads  to 
serious  physical  contradictions.  E.G. ,  we  now  understand  that  (i)  the  clock  retarda¬ 
tion  paradox  would  imply  the  existence  of  absolute  velocities,  which  contradict  the 
postulate  of  the  relativity  of  the  velocities  of  moving  bodies,  and  (ii)  the  assumed 
principle  of  the  invariance  of  the  light  velocity  can  only  refer  (if  at  all)  to  the 
average  light  velocity,  of  a  go-and-return  path  [Builder,  1958a,  b].  Accordingly, 
Einstein's  clock  synchronization  represents  a  “thought-ritual,"  which  has  no  empiri¬ 
cal  value  for  the  measurement  of  the  actual  velocity  of  light  signals  [Alfvdn,  1977], 
and  the  special  theory  of  relativity  is  a  tautology  based  on  average  two-way  signal 
velocities  [Builder,  1958a, b;  Janossy,  1953,  1963;  Ives,  1948]. 

Comprehensive  discussions  of  electromagnetic  substratum  physics ,  from  the 
theoretical  and  experimental  points  of  view,  are  due  to  Janossy  [1953,  1963].  Dirac 
[1958]  and  Kaempffer  [1953]  introduced  the  ether  into  quantum  mechanics.  The 
non-Lorentz  covariant  theory  of  Wilson  [1974]  interprets  elementary  particles  as 
phase  changes  of  an  ether  model.  Winterberg  [1984]  developed  a  nonlinear  relativity 
theory  with  ether  and  a  minimum  length,  which  removes  the  singularities  of  quantum 
electrodynamics  on  a  physical  basis.  It  is  equally  remarkable  that  the  apparent 
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velocity  dependence  of  mass  is  explainable  as  a  particle  interaction  with  the  ether 
[Bagge,  1979]. 

The  Lorentz  transformations  (which  would  hold  in  absence  of  the  ethe'r,  and 
follow  from  the  alleged  invariance  of  the  velocity  of  light  signals),  are  not  the 


sole  transformations  which  leave  Maxwell's  equations  covariant.  Other  (real) 
covariant  transformations  of  Maxwell's  equations  have  been  found  by  Cunningham 
[1909],  Bateman  [1910],  Fushchich  [1978],  and  Fushchich  and  Nikitin  [1982].  Under 
consideration  of  the  quantum-mechanical  commutation  relations,  Winterberg  [1984] 


derived  nonlinear  space-time  transformation  for  high  energy  systems.  Typical  for 
these  transformations  is  the  use  of  conventional  constitutive  relations  D  =  e0E  and 
B  -  u eH  for  the  free  space,  l.e.,  the  vacuum  is  assumed  to  have  invariant  polariza¬ 
tion  properties  [Einstein,  1916,  1921]. 

The  simultaneous  covariance  of  Maxwell's  equations  with  respect  to  Lorentz  and 
other  space-time  transformations  suggests  the  existence  of  covariance  under  an  over¬ 
group  of  these  transformations  [Post,  1962,  1967,  1972,  1979].  The  earliest  predic¬ 
tions  of  the  possibility  of  simultaneous  Lorentz-Galilei  covariance  of  Maxwell's 
equations  go  back  to  Kottler  [ 1922a, b],  Cartan  [1942],  and  van  Dantzig  [1934],  The 
over-group  is  identified  as  Che  nonlinear  set  of  coordinate  substitutions  in  space- 
time  [Post,  1972,  1978].  The  "exclusiveness”  of  Lorentz  covariance  in  relativity 
theory  is  the  artificial  result  of  an  (arbitrary)  restriction  to  the  unimodular 
choice  /-g  *  1  [Einstein,  1916],  which  eliminates  other  transformation  groups  [Post, 
1972,  1978]. 

Kottler  [1922b],  Cartan  [1924],  van  Dantzig  [1934],  and  Post  [1962]  have  shown 
that  the  electromagnetic  field  equations  can  be  brought  into  a  metric-free  represen¬ 
tation,  l.e.,  Maxwell's  equations  exhibit  a  metric-independent  covariance.  Thus, 
Maxwell's  equations  permit  a  manifold  of  space-time  transformations,  if  they  are  not 
restricted  through  the  usual  linear  conaltlutlve  relations  D  *  cE  and  B  *  pH  [Post, 
1972,  1979].  Maxwell's  equations  are  covariant  even  against  transformations  to 

nonlnertlal  frames,  e.g.  rotating  frames  [Post,  1967],  for  certain  constitutive 
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relations  1)  =  F(E)  and  B  =  G(H).  The  accepted  linear  constitutive  equations  appear 
to  be  strictly  valid  only  for  transformations  between  inertial  frames  [Post*  1972, 
1979]. 

Maxwell's  equations  refer  to  an  inertial  frame  V0  in  which  the  electromagnetic 
substratum  is  at  rest,  "w0  =  0,  and  are,  for  this  reason,  not  Galilei  covariant 
[Whittaker,  1951].  We  derive  herein  electromagnetic  field  equations  for  inertial 
frames  V  with  substratum  velocity  w  from  established  axioms  of  physics  (Maxwell  equa¬ 
tions  and  Galileian  relativity  principle).  We  demonstrate  that  these  generalized 
Maxwell  equations,  which  contain  explicitly  the  substratum  velocity  w,  are  covariant 
against  Galilei  transformations .  Since  these  are  transformations  between  inertial 
frames,  the  usual  linear  constitutive  relations  are  assumed  [Post,  1972,  1978]. 

The  electromagnetic  field  equations  with  substratum  flow  vT  and  their  Galilei 
transformations  represent  a  field  theory  for  absolute  or  Galileian  7  ■  (x,y,z)  and  t 
coordinates*  In  inertial  frames  with  substratum  drift  "w,  the  Coulomb  fields  of 
charged  particles  are  deformed  by  convection  so  that  a  physical  length  contraction  of 
material  bodies  in  the  sense  of  Lorentz  [1909]  occurs.  This  length  contraction 
brings  about  a  time  dilatation  since  a  clock  can  be  visualized  as  a  light  signal 
which  is  reflected  anisotropically  in  the  substratum  space  between  two  mirrors  held 
apart  by  a  rod  [Builder,  1958a,  b].  For  this  reason,  also  the  Interrelation  between 
the  absolute  Galileian  and  the  actually  measured  coordinate  and  time  differences  will 
be  discussed,  which  is  determined  by  the  Lorentz  scaling  measure  y  ■  (1  - 

w2/c2)1/2,  i.e.,  a  constant  parameter  for  any  given  inertial  frame  (w). 

The  derivation  of  the  Galilei  covariant  Maxwell  equations  for  inertial  frames 
with  substratum  flow  are  of  interest  for  (i)  mathematical  and  (ii)  physical  reasons. 
The  generalized  electromagnetic  field  equations  provide  physical  foundations  for 
investigations  of  the  electromagnetic  substratum.  As  applications  of  this  theory, 
fundamental  initial-boundary-value  problems  are  solved  analytically  for  electro¬ 
magnetic  signal  propagation  and  induction  in  the  substratum.  We  present  these 
results  for  (i)  theoretical  discussion  and  (ii)  comparison  with  experiments. 
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ELECTROMAGNETIC  FIELD  EQUATIONS  WITH  SUBSTRATUM 

According  to  Maxwell,  Heaviside,  Lorentz,  and  Poincare,  the  electromagnetic 
fields  E°(r°,t0)  and  H°(r°,t°)  in  an  inertial  frame  of  reference  £° (r°, t°,w°  =  0) , 
in  which  the  electromagnetic  ether  is  at  rest,  w°  =  3,  are  determined  by  the 
classical  Maxwell  equations  for  conducting  media  with  velocity  field  [Stratton,  1941]: 


=  -a1°/3t°  ,  (1) 

V°x1l0  =  +3^°/3t°  +  j°  ,  (2) 

V°'^°  =  p°  ,  (3) 

V°-$°  -  0  ,  (4) 

where 

V  -  p°y8  +  a°(t°  +  vV&°)  ,  (5) 

t>*  -  e°t°  ,  V  =■  .  (6) 


As  usual,  the  dielectric  permittivity  ef,  the  magnetic  permeability  y°,  and  the  electric 
conductivity  a0  are  treated  as  isotropic  and  homogeneous.  Equations  (1)  -  (6)  hold  not 
only  for  conducting  media,  but  also  for  insulator,  g  media  including  the  so-called  "va¬ 
cuum''  (  a® s  0)  of  the  system  ,  which  contains  electromagnetic  substratum  at  rest, 
w°  -3.  Equation  (5)  is  Ohm's  law  for  the  current  density  3°  with  space  charge  flowp°v°. 

Let  Eqs.  (1)  -  (6)  be  subject  to  a  Galilei  transformation  of  the  space 
and  time  coordinates  from  the  system  (r° ,  t® ,w°*t))  to  a  system  £(r,t,w  A), 
which  moves  with  the  constant,  but  otherwise  arbitrary,  velocity  u  relative  to 
£®  (0*  of  £*  and  0  of  £  are  assumed  to  coincide  for  t®  *  t  ■  0,  Fig.  1): 

r  -  r°  -  ut°  ,  t  -  t®  ,  (7) 

3 /3 1°  -  3/3t  -  u-7  ,  7®  »  7  ,  (8) 


j  Jt  ~v  x.  <m'j L . A -a  x Mjk' 


with 


F.°  =  E  -  u*  B  , 


-*o  ■+.-*■ 

J  =  j  +  pu 


V  =  V  +  u 


H  =  H 


P°  =  P 


as  Galilei  field  transformations.  Eqs.  (8)  follow  from  Eq.  (7)  by  partial 
differentiation,  3f°  (r°,t°)/3  t®  =  [3  f  (r ,  t) /3  t  ]3 1/3 1®  +  [3 f (r , t) /3r] *3r/3 t° ,  and 
3f®(r®,t®)/3r°  *  [3f(r,t)/3r]‘3r/3r®,  where  3r/3t®  =  -u  and  3r/3r®  =  3  (6._  =  1, 

i  -  J;  -  o,  i  *  j). 

The  Galilei  field  transformations  (9)  -  (11)  are  established  empirically, 
but  will  be  justified  theoretically  by  covariance  requirements.  Furthermore,  since 


medium  density  and  temperature  are  Galilei  invariants. 


e  ■  e 


p  *  w  , 


.  (12) 


The  invariance  of  e  and  p  implies  the  Galilei  invariance  of  the  characteristic 


phase  speed  of  light  [Stratton,  1941] 


o  /  o  o\"l/2  /  ,~l/2 

c  -  (ye*)  «  (ye)  -  c 


Since  the  ether  is  at  rest,  w®  ■  (5,  in  the  system  £®,  the  ether  moves 
with  the  velocity  w  ■  -u  +  ($  in  the  system  £  (Fig.  1).  Application  of  the 
Galilei  transformations  (7)  -  (12)  to  Eqs.  (1)  -  (6)  yields,  therefore,  the 
electromagnetic  field  equations  in  the  inertial  frame  of  reference  £(r,t,w #(5) , 
in  which  the  ether  streams  with  the  velocity  w: 


Vx(!  +  wx|)  -  — (3 / 3 1  +  w*V)5 

7x5  »  +(3/3t  +  w*7)(5  +  ewxf)  +  T  -  Pw 


I, 


/ 


V* (D  +  ew*B)  =  p 

7  •  B  =  0 

where 

-►  -►  -*■  -*  «■*» 
j  »  pv  +  a(E  +  v*B) 

-*■  -► 

D  »  eE  ,  B  =  uH 


(16J) 

(17) 


(18) 

(19) 


Below,  it  will  be  demonstrated  that  Eqs.  (14)  -  (19)  are  covariant  against 
the  three-fold  infinite  number  of  possible  Galilei  transformations.  For  this 
reason,  Eqs.  (14)  -  (19)  are  fundamental  electromagnetic  field  equations,  which 
hold  for  all  inertial  frames  of  reference  £(r,t,w),  in  which  the  ether  flows 
with  constant,  but  otherwise  arbitrary  velocity  w  ^  In  addition,  Eqs.  (14)  - 
(19)  are  approximately  correct  for  inertial  frames  of  reference  £(r,t,w(r,t))  in 
which  the  ether  flow  field  w(r,t)  is  inhomogeneous,  as  long  as  the  spatial  and 
temporal  nonuniformities  of  w(r,t)  have  characteristic  extensions  |Ar|  in  space 
and  j  At |  in  time,  which  are  large  compared  with  the  dimension  and  duration 
of  the  electromagnetic  process,  respectively.  The  latter  conditions  are  frequen¬ 
tly  satisfied  since  most  experiments  are  restricted  to  spatial  regions  which  are 
neglegibly  small  compared  with  the  universe. 

At  interfaces  between  different  media,  boundary  conditions  for  the  electro¬ 
magnetic  field  vectors  ate  required.  These  are  obtained  by  integrating  Eqs. 

(14)  -  (17)  across  the  interface  with  normal  vector  n.  If  n  points  from  medium 
"1"  to  medium  "2"  and  [F]  3  F2  ~  F^»  the  boundary  conditions  in  presence  of 
ether  flow  w  and  interface  motion  v  are: 

n*[E]  -  (£.v)[$]  ,  (20) 

n*[5  +  wx(i$  +  ew*B)  J  «  -(n*v)(5  -f  ewxB]  +  j* 


-  •"  .v.'  A.V..V 


(21) 


s 


n-  [5  +  £w<B]  =  *  ,  (22) 

n-[B]  =  0  (23) 

where  is  the  surface  current  density  and  p*  is  the  surface  charge  density. 

In  the  derivation  of  Eqs.  (20)  and  (21),  the  relations  w*7B  =  -?x(w*B)  and 

-*»  -v  -*■  -v 

wV[D  +  ew*B)  -  pw  *  -7x{w*( D  4-  ew*B)K  both  for  constant  w,  were  used.  The  inter 

action  of  the  ether  with  the  interface  is  negligible,  hence  [w]  =  0.  In  general, 

the  interfaces  separate  regions  of  space  in  which  both  t ,  u,  and  j  are  different, 

and  the  interface  moves  with  a  velocity  field  v.(r,t)  ^  0. 

The  derivation  of  the  generalized  electromagnetic  field  equations  (14)  -  (17 

has  been  carried  through  in  the  MKS  system  for  physical  reasons.  In  comparison 

to  the  cgs  system  (E  and  H  same  units),  the  Giorgian  system  is  superior  since 

(i)  it  introduces  charge  as  a  separate  indeoendent  unit  and  (ii)  it  permits  to 

treat  charge  as  an  invariant  not  only  for  Lorentz  but  for  all  Kottler-Cartan- 

Dantzig-Post  type  space-time  transformations  [Post  1972,  1979]  including  the 

Galilei  transformation.  The  invariance  of  charge  "e"  is  related  to  the  invar i- 
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ance  of  the  characteristic  speed  of  light  c  *  (y  e  )  since  c  -  m  e /e  h 

o  o  o  o  o 

for  dimensional  reasons  (no  mass  dependence,  m^  -  1) .  This  fundamental  relation 
appears  to  indicate  that  electromagnetic  wave  propagation  involves  quantum- 
mechanical  interactions  (h  *  Planck  constant)  of  charges  e  associated  with  the 
ether  medium  (c ,u  ). 


GALILEI  COVARIANCE  OF  ELECTROMAGNETIC  FIELD  EQUATIONS 


The  laws  of  nature  are  of  the  same  form  in  ail  inertial  frames  [Einstein,  1916] 
In  view  of  the  derivation  of  due  electromagnetic  field  equations  (14)  -  (19)  with 
ether  f low  w  they  should  be  covariant  against  Galilei  transformations.  The  rela¬ 
tions  for  the  Galilei  transformation  of  the  independent  variables  and  the 
dependent  fields  from  the  inertial  frame  £(r,t,w)  to  the  inertial  frame  £' (r’.tW)  , 
which  moves  with  the  constant,  but  otherwise  arbitrary  velocity  u  relative  to  £, 
are  given  by  (0  and  O'  of  £  and  coincide  for  t  *  t'  ■*  0,  Fig.  2): 


and 


r*  *  r  -  ut  , 

3/3t  -  3/3t '  -  u-v'  , 

1  «  1'  -  U*B' 


t'  »  t 

7*7' 

H  -  H’ 


I  -  V  +  p'“ 


P  -  p’ 


+  UL 


where 


e  -  e' 


* 


w  *  w  +  u 


■  u  ' 


a  »  o' 


,  (24) 

,  (25) 

,  (26) 

,  (27) 

,  (28) 

(29) 


Substitution  of  Eqs.  (24)  -  (29)  into  Eqs.  (14)  -  (19)  results  in  the  electro¬ 
magnetic  field  equations  for  the  inertial  frame  £'(r',  t',  w')  with  ether  flow  w': 


V'x(1*  +  w'x’S’)  -  -(3/3t'  +  w,-V,)1'  ,  (30) 

V'xfi'  -  +(3/3t’  +  w'-V’MD'  +  E'w'xB')  +  j'  -  p'w’  ,  (31) 

V ' •  (t> '  +  e'w'x'g’)  -  p*  ,  (32) 

V' •!'  «  0  • 


where 


p' v'  +  o ' (t '  +  v'xS') 


V 


(34) 


ft'  *  e'E'  ,  B'  =  u'H* 


(35) 


The  Eqs.  (30)  -  (35)  for  the  inertial  frame  ^'(r'.t'.w*)  are  indeed  of  the  same 
form  as  the  Eqs.  (14)  -  (19)  for  the  inertial  frame  £(r,t,w).  Thus,  the  Galilei 
covariance  of  Che  proposed  electromagnetic  field  equations  (14)  -  (19)  with 
ether  flow  w  is  demonstrated.  This  covariance  is  the  necessary  condition  for 
Eqs.  (14)  -  (19)  to  be  generally  valid. 

It  is  remarkable  that  the  Galilei  transformations  (26)  of  the  fields  E 
and  ft  are  not  symmetric.  In  view  of  Eq.  (29),  Eq.  (22)  can  also  be  stated  as 

ft  «  ft’-  e'u* B'  ,  ft  -  B'  .  (36) 

The  analogous  transformation  formula  for  the  magnetic  field,  B  -  B'+  u’ftxft',  does 
not  render  Eqs.  (14)  -  (19)  Galilei  covariant.  This  relation  is  incompatible 
because  it  implies  magnetic  charges,  V*B  4  0,  for  which  no  experimental  evi¬ 
dence  exists. 


According  to  Eq.  (16),  an  ether  flow  w  transverse  to  the  magnetic  field 
Induces  an  electric  charge  density  Pw, 

VD  -  p  +  p  ,  P  -  -V*  (ew*ft)  .  (37) 

w  w 

This  interesting  effect  could,  in  principle,  be  used  to  detect  the  ether  flow  w 
through  space  charge  measurements.  Unfortunately,  it  appears  that 

P  -  c_2w*V*H  (38) 

w 

is  very  small  for  laboratory  experiments  (w  <  10^m/sec,  V*H  <  lO^A/m^)  and 
even  for  cosmic  situations,  e.g.,  for  quasars  (w  %  c,  ?*ft  <  10  ^A/m^). 


w+ 


WAVE  EQUATIONS  FOR  ELECTROMAGNETIC  POTENTIALS 

The  basic  Eqs.  (14)  -  (19)  represent  a  system  of  coupled  partial  differ¬ 
ential  equations  of  first  order  for  the  electromagnetic  field  E(T,t)  and 
1i(r,t).  With  regard  to  mathematical  applications,  it  is  desirable  to  work  with 
decoupled  wave  equations  for  the  scalar  potential  $(r,t)  and  the  vector  poten¬ 
tial  ^(r,t).  From  these,  the  electromagnetic  fields  are  derived  as  partial 
derivatives, 

1  -  -V*  -  al/3t  ,  8  =  VxX  .  (39) 

Since  an  arbitrary  vector  field  A(r,t)  consists  of  a  solenoidal  and  an  irrota- 
tional  part,  and  V*A  is  still  unspecified  by  Eq.  (39),  a  is  subject  to  the 
"ether  gauge" 

V*^  *  -ye  (3/3 1  +  w*V)($  -  w*"&)  ,  (40) 

which  reduces  to  the  Lorentz  gauge  for  w  *  t).  The  relations  in  Eq.  (39) 
satisfy  Eqs.  (14)  and  (17)  identically.  Elimination  of  1)  ■  and  ^  »  pfi 
from  Eqs.  (15)  and  (16)  by  means  of  Eq.  (39)  yields,  under  consideration  of 
Eq.  (40),  the  inhomogeneous  wave  equations  for  the  vector  potential  X(r,t)  and 
the  scalar  potential  y(r,t)  •  $(r,t)  -  w*^(r,t): 

[ye(3/3t  +  w*7)2  -  72$  *  y  (}  -  pw)  ,  (41) 

[ye(3/3t  +  w*7)2  -  72]($  -  w*X)  *  p/e  ,  (42) 

_2 

where  ye  ■  c  .  These  hyperbolic  equations  exhibit  the  convective  influence 
of  the  ether  flow  w  on  the  electromagnetic  potentials.  Eqs.  (41)  and  (42) 
reduce  to  the  conventional  wave  equations  for  ^(r,t)  and  $(r,t)  in  the 
special  case  of  the  Maxwell  frame  with  resting  ether,  w  •  t)  [Stratton,  1941]. 


the  initial  and  boundary  conditions  for  A(r,t)  and  7(r,t)  =  i> <r , t )  - 
w*i£(r,t)  follow  from  those  for  E(r,t)  and  B(r,t)  via  Eq.  (39).  Thus,  the 
solutions  ^(r,t)  and  ’?(r,t)  of  the  wave  equations  (41)  and  (42)  can  be  estab¬ 
lished  for  known  current  and  space  charge  distributions  l(r,t)  and  p(r,t).  As 
an  example,  we  give  the  solutions  of  fundamental  retarded  potentials  of 
Eqs.  (41)  and  (42),  which  are  generated  by  the  sources  l(r,t)  and  p(r,t): 


^(r,t)  =  j-  /If R_1[j (r*,t-R/c)  -  wp (r*,t-R/c) ]d^r* 


2 

$(r,t)  *  t~  /// R  1  [  (1  -  ^r)p  (r*,t-R/c)  +  c“2w*  j  (r* ,  t-R/c)  Jd3?*  (44) 


where 


/(r  -  r*)2  ,  R(r,r*)  -  R!(r',r*') 


,  (45) 


is  the  distance  between  the  field  point  (r)  and  the  source  point  (r).  By 
Eqs.  (43)  and  (44),  the  sources  j(r  ,t)  and  p (r  ,x)  at  the  source  point  r  con¬ 
tribute,  at  the  retarded  time  x  ■  t  -  R/c,  to  the  fields  X(r,t)  and  $(r,t)  at 
the  point  r  and  time  t  of  observation.  Since  R(r,r*)  and  c  *  (ue)  are 
Galilei  invariants  by  Eqs.  (24)  and  (29),  the  time  retardation  At  ■  -  R/c  is  a 
Galilei  invariant  as  expected  [At  *  At*  by  Eq.  (24)]. 

The  retarded  potential  solutions  (43)  and  (44)  refer  to  the  inertial 
frame  £(r,t,w),  in  which  the  ether  velocity  is  w  j*  They  indicate  that  the 
corresponding  solutions  for  the  Maxwell  frame  £°(r#,t#,w#  *^$)  are  approximately 
valid  also  for  |w®|  <<  c.  However,  for  |w|  -*•  c,  the  effects  of  the  ether 
flow  w  on  the  potentials  A(r,t)  and  $(r,t)  ace  quantitatively  significant. 

Note  that  also  the  nature  of  the  $(r,t)  solution  changes  with  respect  to  its 
sources  as  |w|  increases  from  0  towards  c. 


1  r-rww v nfii ,■  <a  ■» w i mb w, uum.  v f  jirji va^ne 


13 


The  ether  gauge  (40)  and  the  wave  equations  (41)  and  (42)  are  covariant 
against  the  Galilei  transformations  (24)  -  (29).  This  covariance  is  obvious. from 
the  co-/  invariance  of  the  operators  and  field  expressions  in  Eqs.  (40)  -  (42), 


namely : 

3/3t  +  w*7  =  3/3t’  +  w '  •  7 '  ,  7  -  7’  ,  (46) 

t=P  ,  (47) 

<t  -  W’X  *  $'  -  w' •%'  ,  (48) 

J  -  pw  *  j *  -  p 'w'  ,  (49) 

P  -  P '  ,  (50) 


and  Eq.  (29).  The  unpriraed  and  primed  fields  refer  to  the  inertial  frames 
£(r,t,w)  and  J’(r ' ,  t ' ,  w 1 )  ,  respectively.  Since  the  relative  system  velocity  is 
£  »  w  -  w'  *  v  -  v'  (Fig.  2),  Eqs.  (46),  (48),  and  (49)  imply  the  Galilei  co-/ 
invariants 

3/3t  +  v*7  -  3/3t' +  v'*7'  ,  (51) 

♦  -  v-A  ■  -  v ' • A'  ,  (52) 


j  -  pv  ■  j*  -  p  'v' 


(53) 


In  terms  of  the  electromagnetic  fields,  the  Galilei  invariants  (47),  (48),  and 
(52)  assume  by  Eq.  (39)  the  form  [Stratton,  1941]. 


1  -  S’ 


1  +  ■  "S*  +  w'x^* 


^  +  vx?  ■  +  v’xt’ 


(54) 

(55) 

(56) 
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By  means  of  the  invariants  (46),  (49),  (50),  and  (53)  -  (56),  the  Galilei 

covariance  of  the  basic  electromagnetic  field  equations  (14)  -  (19)  with 

ether  flow  w  is  now  immediately  recognized. 

In  the  special  case  w  »  (5,  the  new  wave  equations  (41)  and  (42)  with 

ether  flow  combine  to  the  relativistic  wave  equation  t3{A,i$/c}  *  -u{j,icp} 

2  -2  o 

with  the  Lorentz  covariant  operator  □  ■  V  -  c  S“.  This  D'Alembertian  can, 
in  no  way  whatever,  describe  anisotropic  light  propagation  or  the  nonreci¬ 
procal  asymmetry  between  the  clockwise  and  counterclockwise  beams  observed 
in  the  Sagnac  experiment  [Post,  1967].  Anisotropy  and  nonreciprocity  require 
mixed  space-time  derivatives  7^  in  the  wave  equation  [Post,  1967],  These 
are  generated  by  the  ether  convection  w*V  in  the  proposed  wave  equations  (41) 
and  (42)  with  the  nonsymmetric  operator  □  »  7  -  c  (3t  +  w*7)  .  The  space- 
time  symmetry  is  destroyed  by  the  ether  flow  w,  i.e.,  exists  only  in  the 
ether  rest  frame  where  w  a  ft. 


ELECTRODYNAMIC  CONSERVATION  THEOREMS 


In  the  generalized  electrodynamics  with  ether  flow  w,  charge,  energy, 
and  momentum  of  the  fields  are  conserved  in  all  inertial  frames  £(r,t,w)  with 
ether  motion  w.  These  conservation  theorems  are  presented  and  shown  to  be 
Galilei  covariant. 

Charge  Conservation.  The  divergence  of  Eq.  (15)  yields,  under  consideration 
of  the  vector  identity  for  arbitrary  fields  a(r,t), 

V*  [(3/3t  +  w-V)a]  -  (3/3t  +  w*V)7-a  ,  (57) 

and  Eq.  (16),  the  conservation  equation  for  the  charge  density  in  the  inertial 
frame  £(r,t,w): 

(3/3t  +  w-7)p  -  -7-tf  -  pw)  .  (58) 

This  equation  is  Galilei  covariant  in  view  of  Eqs.  (46),  (49),  and  (50). 

Eq.  (58)  is  equivalent  to  the  usual  charge  continuity  equation  [Stratton,  1941 J 

3p/3t  -  >7*1  *  .  (59) 

Electromagnetic  Energy  Conservation.  In  accordance  with  the  vector  relations 
for  arbitrary  fields  a(r,t)  and  ft(r,t), 

V*(axt)  -  $  •Vxa  -  a*Vxt  ,  (60) 

a*(3/3t  +  w*7)a  -  (3/3 t  +  w-7)a2/2  ,  (61) 

scalar  multiplications  of  (i)  Eq,  (14)  by  ft  and  (ii)  Eq.  (15)  by  (1  +  wxft) , 

-  and  subtraction  of  the  resulting  equations  (i)  and  (il)  -  ,  results  in  the 
conservation  equation  for  the  field  energy  in  the  inertial  frame  £(r,t,w): 
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(3  /3  t  +  w*7)  [y(E  +  wx  B)2  +  y  H“]  =  -7 • [ (E  +  wxH)x"fi]  -  (j  -  pw)  •  (£  +  wxB)  . 

(62) 

This  equation  is  Galilei  covariant  by  Eqs.  (46),  (49),  (54),  and  (55). 


Equation  (62)  can  be  rewritten  in  the  form 

(3/3t  +  w-V)(U  +  U  )  =  -7-1  -  Q  (63) 

e  m 

where 

Ug  «  ye  (t  +  wxt)2  ,  (64) 

um  *  y u  tl2  ,  (65) 

^  -  ( E  +  wx  B)x H  ,  (66) 

Q  “  (1  -  pw).(t  +  w*l)  ,  (67) 


are  the  electric-  and  magnetic  energy  densities,  the  Poynting  vector,  and  the 
ohmic  power  density  in  presence  of  ether  flow  w.  Note  that  Eqs.  (64)  -  (67) 
define  generalized  concepts  for  energy  density,  energy  flux,  and  power  density 
of  the  electromagnetic  field,  which  are  Galilei  invariant. 
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The  dyadics  (71)  and  (72)  are  Galilei  invariants  by  Eqs.  (54)  and  (55). 
Vectorial  multiplications  of  (i)  Eq.  (14)  by  *c(E  +  w*B)  and  (ii)  Eq.  (15) 
by  and  addition  of  the  resulting  equations  (i)  and  (ii)  ,  gives  the  conser- 

r*  — ►  — ► 

vation  equation  for  the  electromagnetic  momentum  in  the  inertial  frame  >(r,t,w): 


(3/3 1  +  w»V) {ye (E  +  wxI)xH]  =  V  •  +  T )  -  o  (E  +  wxB)  -  (j  -  pw)«B  ,  (73) 

e  m 


(3/3 1  +  w-V)(ue?)  =  +  T  jxB 

e  m 


•  (74) 


These  equations  are  Galilei  covariant  by  Eqs.  (46),  (49),  (50),  (55),  (71), 
and  (72).  Equation  (73)  or  (74)  is  a  relation  through  which  the  forces  on  charge 
p(r,t)  and  current  j(r,t)  densities  can  be  expressed  by  the  electromagnetic  fields 
^(r,t)  and  B(r,t)  in  the  medium,  which  is  assumed  to  be  homogeneous  with  respect 


to  e,  u,  and  a.  In  particular,  if  the  electromagnetic  momentum  change  (3^  +  w?) 

(pc?)  is  negligible,  the  electromagnetic  force  density  equals  the  divergence  of 

the  field  stress  tensor,  oE  +  jxB  »  +*?  )  [Stratton,  1941]. 

e  m 

Within  the  frame  of  the  Galilei  covariant  electrodynamics  for  inertial 
systems  £(r,t,w)  with  ether  flow  w,  generalized  conservation  theorems  for  energy 
[Eq.  (67)]  and  momentum  [Eq.  (73)]  of  the  electromagnetic  field  have  been  found. 
These  reduce  exactly  and  approximately  to  the  ordinary  electromagnetic  conser¬ 
vation  equations  for  ether  velocities  w  -  t)  and  |w|  <<  c,  respectively.  Equa¬ 
tions  (67)  and  (73)  predict  significant  physical  effects  for  large  ether  veloci¬ 
ties  |w|  ■+  c,  e.g.,  in  the  vicinity  of  distant  galaxies  and  quasars  with  extremely 
large  ether  expansion  drifts.  A  verification  of  the  ether  terms  in  Eqs.  (67) 


and  (73)  by  means  of  laboratory  experiments  may  be  difficult,  since  macroscopic 

A 

platforms  can  presently  not  be  accelerated  to  speeds  v  >  10  m/sec. 
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SIGNAL  PROPAGATION  IN  SUBSTRATUM 

According  to  the  solutions  of  Maxwell's  equations  (1)  -  (6)  for  the 

inertial  frame  of  the  resting  ether  (w°  =  t)) ,  light  signals  propagate  with  an 

-1/2 

isotropic  speed  v®  ■  c  a  (pc)  ,  which  is  independent  of  the  velocity  and 
acceleration  of  the  source.  .  This  independence  from  velocity  and  acceleration 
of  the  emitter  holds  also  for  signal  propagation  described  by  the  Galilei  co¬ 
variant  electromagnetic  field  equations  (14)  -  (19).  However,  the  signal  velocity 

v  =■  v  (r,w)  is  no  longer  isotropic  in  inertial  frames  Y(r,t,w)  with  ether  flow 
s  s 

w  +  t). 

In  order  to  understand  signal  propagation  in  the  ether,  consider  the 
elementary  excitation  of  an  electromagnetic  wave  pulse  by  the  sudden  application 
of  a  current  pulse' 3*(t)  [A/m]  to  the  surface  x  -  0  of  an  ideal  conductor  (a  ■+•  ®) 
at  time  t  *  0, 

j*(t)  -  (0,0,JH(t)}  ,  x-0  ,  t>0  ,  (75) 

where  H(t)  is  the  Heaviside  step  distribution,  dH(t)/dt  *  6 (t) .  The  resulting 
electromagnetic  wave  emitted  from  the  "sheet  antenna"  at  x  ■  0  at  time  t  *  0  is 
of  the  form  (Fig.  3) 

t  -  {0,B(x,t),0}  ,  1  -  (0,0,E(x,t)}  .  (76) 

The  propagation  of  the  wave  (76)  in  the  charge  (p  »  0)  and  current  (3  -  0)  free 
space  x  >  0  with  ether  flow  parallel  to  the  wave  propagation  (Fig.  3) 

w  -  {w(1,0,Q}  ,  w()  i  0  ,  (77) 


19 


(3  /3 1  +  W||  3  /3  x)  2B  =  c232B/3x2 

9 

0  < 

x  <  =° 

9 

(78) 

B(x,t*0)  =  0 

9 

0  < 

X  <  ® 

9 

• 

(79) 

3B(x,t=0)/3t  ■  0 

9 

0  < 

X  <  ® 

9 

(80) 

B(x=0, t)  -  u(l  -  «j2/c2)_1j*(t) 

9 

0  < 

t  <  00 

« 

(81) 

The  Laplace  transform  method  gives  the  solution  of  Eqs.  (78)  -  (81)  as  the  integral 
functional 


where 


a+i®  ,  .  , 

c,t)  *  977  /  F(s)ets-<*-”lltWcds 


f  F(s)e(14^l  /C)t3ds  -  p(l  -  «)f/c2)‘1j*(t) 


by  the  boundary  condition  (81) .  The  integral  equation  (83)  for  the  Laplace 
amplitude  F(s)  has  closed-form  solutions  for  certain  j*(t)  functions,  e.g. , 


F (s)  -  u(l  -  w,2/c2)  ha  1  , 


w„  <  c 


,  (84) 


for  the  boundary  value  (75).  Equations  (82)  and  (84)  combine  to  the  wave  pulse 
solution : 


B(x,t)  -  p (1  -  w2/c2)’13H(t  - 


,  x  >  0 


E(x,t)  -  -(^|  +c)u(l  -  w2/c2)  ^(t  -  ,  x  >  0  ,  (8 

by  Eq.  (14).  Equations  (85)  and  (86)  follow  from  the  discontinuous  integral 
[Eqs.  (82)  and  (84)], 


—  /’0l+i*g-1ets-<x-w||t)s^c 

J  9  * 


x-w,,  t 

ds  «  H(t - — — ) 


-  •  *.**V  V«'\>  *>  *  jfc 


where 


x-w„  c 


H(t - —-)  =  H(t 


w„  +c 


<  c 


(88) 


The  results  (85)  and  (86)  are  physically  interesting  in  several  respects. 

As  typical  for  hyperbolic  equations  or  processes  with  finite  characteristic  speed 
c  <  °>,  the  wave  fields  are  discontinuous  at  the  wave  front  (Fig.  3).  The  posi¬ 
tion  and  speed  of  the  signal  front  are  in  the  inertial  frame  £(x,y,z,t,w„  ) 


5(t)  =*  (w||  +  c)t  ^  ct  ,  w|r  ^  0 

i  a  0 


(89) 

(90) 


Accordingly,  the  electromagnetic  signal  propagates  with  a  velocity  v*  i  c  depend- 

s 

ing  on  whether  the  propagation  occurs  in  (>)  or  oppo»At  (<)  to  the  direction  of 

the  ether  flow  Wjj  .  Only  in  inertial  frames  moving  with  the  ether,  the  signal 

propagates  with  the  velocity  v#  ■  c  in  all  directions,  i.e.,  the  Galilei  invariant 

s 

-1/2 

c  =  (ye)  is  the  speed  of  light  relative  to  the  ether.  Apparently,  c  repre¬ 
sents  the  upper  limit  for  the  ether  velocity,  |w|  <  c,  since  |b[  »  and  fl|  -*•  ® 

for  |w|  -*•  c  by  Eqs.  (85)  and  (86).  It  is  also_  seen  that  for  the  electric  and 
magnetic  field  intensities  of  wave  phenomena 

|£|  “  lw||  +  CH^I  »  |W||  +  c|  «  c  ,  W||  0  .  (91) 


This  relation  could  possibly  be  used  for  the  determination  of  the  ether  speed 
|W||  |  through  very  accurate  interferometric  comparisons  of  the  electric  and  magne¬ 
tic  amplitudes  of  waves. 

In  the  same  way,  the  initial-boundary-value  problems  for  signals  propagat¬ 
ing  in  directions  perpendicular  to  the  ether  velocity  w  can  be  solved,  with  v^  » 
2  2  1/2 

(c  -  w  )  '.  Accordingly,  electromagnetic  signal  propagation  is  anisotropic 

v*  i*  v1  in  inertial  frames  with  ether  flow. 

s  s 


ELECTROMAGNETIC  INDUCTION  IN  SUBSTRATUM 


As  an  elementary  model  for  the  generation  of  electromagnetic  pulses  by 
induction  in  the  ether,  consider  a  plane  copper  slab  (initially  at  x  =  0) ,  which 
is  accelerated  at  time  t  =  0  across  a  homogeneous  magnetic  field  Bq  =  {0,Bo>0} 
to  a  velocity  v  *  {a(t),0,0},  so  that  the  position  of  the  front  surface  is  at 
x  =  a(t)  at  time  t  with  a(t=0)  =  0  (Fig.  4).  The  duration  At  of  the  (explosion 
driven)  piston  motion  is  assumed  to  be  small  compared  with  the  field  diffusion 
time  tD  »  uad  (slab  thickness  d)  .  The  motion  of  the  quasi-ideal  conductor  (a-*00 
across  induces  at  its  front  surface  an  electromagnetic  field  of  the  form 

15  *  {0,- 3A(x,t)/3x,0>  ,  £  *  {0,0,-3A(x,t)/3t}  (92) 

which  propagates  into  the  space  x  >  a(t)  of  Fig.  4.  Let  the  ether  flow  be 
parallel  to  the  direction  of  wave  propagation,  w  *  ,0,0},  w^  ^  0.  The  vector 

potential  A  *  (0,0,A(x,t)}  is  then  determined  by  the  hyperbolic  initial-boundary- 
value  problem  for  the  homogeneous  wave  equation  (41): 

(3/3 1  +  wj|3/3x)2A  »  c232A/3x2  ,  a(t)  <  x  <  «  ,  (93) 

A(x,t-0)  -  -Bqx  ,  0  <  x  <  ®  ,  (94) 

3A(x,t-0)/3t  -0  ,  0  <  x  <  co  ,  (95) 

[3A(x,t)/3t  +  a(t)3A(x,t)/3x]x=a^t^  -  0  ,  0  <  t  <  »  .  (96) 

Equations  (94)  -  (96)  consider  that  B(x,t-0)  -  Bq,  E(x,t-0)  -  0,  and  E(x,t)  + 
v(t)B(x,t)  *  0  at  the  moving  piston  (a  *  <*>)  surface  x  ■  a(t).  The  solution  of 
Eqs. (93)  -  (96)  is  by  the  Laplace  transform  method 
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ct+i® 


AOc.t)  •  -Box  +  jJj-/  .  C(S)ets'(*^lt)s/<:ds 


Oi-ioo 


where 


A  /“+i"S<:(S,ets-'a<t>-'lltls/'=ds  -  a(t)B  /[l  - 
2ir  i  J  .  o  c 

a-i® 


a(t)  -  w„ 


(97) 

(98) 


by  the  boundary  condition  (96) .  The  integral  equation  (98)  for  the  Laplace 
amplitude  G(s)  has  closed-form  solutions  for  certain  piston  motions  a(t),  e.g. , 
for  the  Heaviside  velocity  pulse 


a(t)  =  vQH(t) 


G(s) 


v  B  s~2/(l  - 
o  o 


<  c 


(99) 


From  Eqs.  (97)  and  (99)  result,  under  consideration  of  Eq.  (92),  the  solution 
for  the  induced  electromagnetic  pulse  fields: 


(l  +  w../c) 

A(x,t)  *  -B  x  +  : - ; - r~\  v  B  (t - — )H(t 

*  O  (1  /C  -  VQ/c)  0  O'  Hjj  +c 


---  --)  ,  X  >  V  t  , 

US.  +  c  o 


B(x,t)  ■  B  + 


v  /c 
o 


+  . - rrr  B  H(t - — -) 

O  (1 /C  -  Vo/c)  O'  Wj|  +C 


’  x  <  vot  . 


,  *  >  v  * 


(100) 


vo/c 


E(x’C)  *  ‘  (1  +  «,  7c  -v  7c)  *1l  +  c>BoH(t  ' 


-  0 


*  *  <  V  ’ 


,  x  >  vot  , 


,  X  <  v  t  . 
o 


(101) 


(102) 


In  the  derivation  of  Eqs.  (100)  -  (102),  the  discontinuous  Laplace  integrals 
have  been  expressed  as  Heaviside  impulse  functions,  in  accordance  with  Eqs.  (87) 
and  (88) . 

Equations  (101)  and  (102)  indicate  that  the  electromagnetic  wave  pulse 
induced  at  the  front  surface  x  -  vQt  of  the  conductor,  which  pushes  the  flux  of 
the  magnetic  field  Bq,  is  discontinuous  at  the  wave  front  S(t)  *  (wy  +  c)t  ^  ct 
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for  wjj  ^  0,  which  propagates  with  the  speed  vg  “  wy  +  c  ^  c  for  wj|  5  0  in  the 
inertial  frame  £  (x,y,z,t,wjj  ) .  Figure  4  shows  the  induced- magnetic  field  pulse 
B(x,t)  *  B(x,t)  -  Bq,  which  occupies  the  space  vQt  <  x  <  (wjj  +  c)t  at  time  t. 
Thus,  the  electromagnetic  signal  propagates  with  a  speed  vj,  ^  c  in  the  inertial 
frame  £  depending  on  whether  the  propagation  occurs  in  (>)  or  opposite  (<)  to 
the  direction  of  ether  flow  w.  The  electric  and  magnetic  wave  intensities  are 
interrelated  by  { E [  *  f +  c|*|b|.  Equations  (101)  and  (102)  indicate  that  for 
ether  drifts  W||  opposite  to  the  piston  velocity  vq 


B(x,t)  -  Bo[  -  - 

for 

vo  "  w||  -  c 

o 

V 

f 

f 

(103) 

E(x,t)  |  -*■  - 

for 

vQ  -  w„  -  c 

•  "il  ' 0 

• 

(104) 

Accordingly,  copious  amounts  of  radiation  would  be  produced  by  magnetic  flux  (Bq) 
pushers  in  regions  of  space  with  large  ether  drifts  10-1c  <  | w^  |  <  c  (distant 
galaxies,  quasars).  On  principle,  even  the  limit  vq  —  wj|  -  c  is  achievable  if 
-W||  deviates  from  c  by  not  more  than  the  conductor  speed  vq.  It  is  remarkable 
that  the  limit  vq  —  Wj|  ■  c  is  not  even  approximately  realizable  in  absence  of 
ether  flow  (w„  »  0)  since  v  «  c  for  macroscopic  bodies. 


•»’  ♦« 


DISPERSION  AND  FREQUENCY  SHIFT  IN  SUBSTRATUM 


In  a  homogeneous  dielectric  (e,y)  or  the  free  ether  space  (£0>y0)  the 


Maxwell  frame  £°(r0,to,w°  =  5)  in  which  the  ether  is  at  rest,  the  frequency  and 
wavelength  of  electromagnetic  waves  are  interrelated  by  w  =  j tc [  c  where  |^|  *  2tt/X, 


In  an  inertial  frame  £(r,t,w)  with  ether  flow  w,  the  dispersion  of  monochromatic 


waves  of  frequency  w 


^(r.t)  *  Re^e 


i(1c*  r-yt) 


(105) 


in  absence  of  space  charges  (p  ■  0)  and  electric  currents  (j  *  t5)  is  determined 


by  the  homogeneous  wave  equation  (41)  for  the  vector  potential  with  ether  flow  w. 


(3/3 1  +  w-V)2^  -  c2?2! 


(106) 


Substitution  of  Eq.  (105)  into  Eq.  (106)  yields,  since  ^  (condition  for 


nontrivial  solution),  the  dispersion  equation  w  =  w(k)  for  electromagnetic  waves 


in  the  ether: 


w  •  |k|c  +k*w  ,  c  *  (ye) 


(107) 


Thus,  in  presence  of  ether  flow  w,  electromagnetic  waves  exhibit  in  the  inertial 


frame  £  a  frequency  shift 


Aw  ■  1t*w  ^  0  ,  cos(k,w)  %  0 


(108) 


which  vanishes  for  propagation  k  perpendicular  to  the  ether  flow  w.  Equation  (108) 


predicts  a  red-shift  (Aw  <  0)  or  a  blue—  shift  (Aw  >  0)  for  electromagnetic  waves 


propagating  with  a  wave  vector  component  kj|  which  is  antiparallel  or  parallel  to 


•  i  < .V*  <«  V.  i*.  V.  W 


Lat:  4.*  |Lr  4  iu**.*4  A. a 


the  ether  velocity  w.  The  physical  reason  for  this  frequency  shift  is  the  propa- 
gation  of  waves  and  signals  with  the  phase  velocity 


v  *  c  +  (£'w)/]1c| 


(109) 


in  the  inertial  frame  £,  i.e.,  the  propagation  of  electromagnetic  waves  with  the 
velocity  c  relative  to  the  moving  ether.  For  laboratory  and  terrestrial  applica¬ 
tions,  the  frequency  shift  is  in  general  small 


Aw  |  «  a*  for  |  w |  «  c 


(110) 


The  red-shifts  of  spectral  lines  observed  in  distant  galaxies  have 
been  the  subject  of  different  explanations,  with  the  Doppler  red-shift  due 


to  the  rapid  recession  of  the  galaxies  [Humason,  1956] 


being  now 


widely  accepted.  •  The  longitudinal  Doppler  red-shift  A(j>d  *  -|k[*|Vp|,  where 
Vq  is  the  velocity  of  the  light  source  (galaxy) ,  is  of  the  same  order-of-magnitude 
as  the  longitudinal  ether  red-shift  from  Eq.  (108),  Au>  *  -|k|*|w|,  if  w  «  Vp. 

Since  the  ether  probably  "expands"  as  the  masses  of  the  universe  recede  ( relative 
to  theearth),  both  the  Doppler  and  ether  red-shifts  have  to  be  considered  in  the 
evaluation  of  the  velocity  v^  »  w  of  the  galaT.ies  from  the  experimental  red-shift 
data30  cAX/X.  These  indicate  ether  drift  velocities  w  -  lO^c  to  4xl0_1c  [Humason,  1956] 
The  ether  red-shift  represents,  therefore,  another  physical  effect  for  the 
experimental  Investigation  of  the  substratum. 

The  presented  applications  provide  an  impression  of  the  electromagnetic 
phenomena  which  can  be  expected  in  presence  of  ether  flow.  The  theory  has  signif¬ 
icant  other  applications  in  connection  with  the  interaction  of  radiation  and 
charged  particles  with  the  ether  under  laboratory  conditions. 


.N  V>- •»  JVkV  .'-iX.-- 


GALILEAN  AND  MEASURED  COORDINATES 

The  generalized  electromagnetic  field  equations  (1A)-(19)  with  substra¬ 

tum  and  their  Galilei  transformation  (24)-(28)  represent  a  field  theory  in  terms 
of  absolute  or  Galilean  space  r  *  (x,y,z)  and  time  t  coordinates,  i.e.,  the  Gal- 
ilean  coordinate  differences  are  the  same  in  all  inertial  frames,  Ar  -  inv  and 
At  *  inv.  According  to  Lorentz  [1909],  a  measuring  rod  resting  parallel  to  the  ether 
velocity  w  in  an  inertial  frame  V(r,t,w)  has  there  the  reduced  length  L(w)  = 

L0(l  -  w2/c2)l/2  due  to  the  flattening  of  its  microscopic  Coulomb  fields 
by  the  ether  flow  (Lq  »  proper  rod  length  in  the  ether  rest  frame 

Recognizing  that  a  clock  can  be  visualized  as  a  system  reflecting  a 
light  signal  back  and  forth  between  two  mirrors  held  apart  by  a  rod.  Builder 
[1958  a,b  ]  demonstrated  that  the  period  of  a  clock  is  increased  to  T(w)  * 

To/d  “  w2/c2)l/2  in  the  Inertial  frame  J(r,t,w)  with  ether  flow  w,  as  the 
combined  result  of. the  rod  contraction  and  the  anisotropy  of  light  propagation 
between  the  mirrors  (T0  »  proper  clock  period  in  the  ether  rest  frame 

For  these  physical  reasons,  the  differences  A  of  the  absolute  Galilean 
coordinates  r,t  of  the  inertial  frame  £(r,t,w)  are  related  to  the  measured  (m) 
space  and  time  coordinate  differences  in  this  reference  frame  by  ( n  and  1  to  w) 

Ar"  -  yAr "  ,  Ar*  -  Ar*  ,  At  -  At/y  ,  (111) 

m  m  m 

where 

y  ■  (1  -  w2/c2)1^2  ,  0  <  Jw|  <  c  .  (112) 

The  measured  coordinates  rjj,  and  tm  vary  in  accordance  with  the  y(w)  of  the 
respective  inertial  frame  £.  For  mathematical  (y-1  <  »)  and  physical  (violation 
of  causality  principle)  reasons,  the  ether  speeds  are  restricted  to  values 
0  <  Jw  j  <  c  (in  agreement  with  observation).  Since  the  scaling  factor  y(w)  is  a 
constant  for  a  given  inertial  frame  £(r,t,w),  the  calculation  of  the  measured 


(m)  coordinates  from  the  Galilean  coordinates,  and  vice  versa,  by  means  of 
Eqs.  (Ill)  -  (112)  is  elementary  in  applications. 

In  order  to  illustrate  the  use  of  the  Galilean  space  and  time  coordi¬ 
nates,  the  Doppler  effect  shall  be  analyzed  (within  the  frame  of  the  presented 
theory)  for  an  emitter  A  and  a  receiver  B,  which  move  along  the  x-axis  with 
the  velocities  u  5  0  and  v  5  0  relative  to  the  ether  frame  £  ,  respectively. 
The  Galilean  position  coordinates  of  A  and  B  are  at  the  absolute  time  t  >  0 

of  the  ether  frame  V 

-o 

xA(t)  -  a  +  ut  ,  xg(t)  =  b  +  vt  ,  (113) 

where  a  and  b  are  the  initial  positions.  The  Galilean  time  t  is  counted  by 
similar  (synchronized)  clocks,  distributed  over  the  ether  frame  (isotropic 
light  propagation) .  The  time  periods  of  the  emitter  and  receiver  shall  be  the 
same,  TA  -  T  *  Tfl,  when  A  and  B  are  at  rest  in  Tq.  Let  A  emit  signals  with 
period  TA  at  the  times 

t  *  t  +  nT.  ,  n  ■  0,1,2,  •  •  •  .  (114) 

no  a 

At  a  time  t  >  t  ,  the  n-th  signal  emitted  by  A  reaches  the  location  on  the 
n 

x-axis 

I 

xn(t)  -  xA(tn)  +  c(t  -  tn)  -  a  -  (c  -  u)tn  +  ct.  (115) 

Accordingly,  the  n-th  signal  of  the  emitter  A  will  hit  the  receiver  B  at  a 

time  t*  determined  by  x  (t*t*)  ■  x_(t*t*),  or 
n  n  n  d  n 

a  -  (c  -  u)t  +  ct*  ■  b  +  vt*  .  (116) 

n  n  n 


Hence 
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where 

t*  s  tQ  +  (b  -  a)/(c  -  v)  , 


.(118): 


T*  a  ta(c  '  u)/(c  '  v)  *  (119) 

Equation  (119)  Indicates  that  the  period  of  the  emitter  A  is  "observed"  as 

a  period  T*  at  the  location  of  the  receiver  B.  So  far,  the  considerations 

A 

have  been  purely  "classical"  or  "Galilean." 

g 

In  experiments,  one  compares  the  period  (measured  at  B)  of  the  receiver 
B  with  the  period  T*^  (measured  at  B,  too)  of  the  emitter  A.  According  to  Eqs. 
(Ill)  and  (112) 

TJm  '  *J/<1  -  “V>1/2  .  V  V(1  '  ''2/=2)1/2  U20) 

since  A  moves  with  the  volocity  u  and  B  moves  with  the  velocity  v  relative  to 
the  ether,  T  .  By  Eqs.  (119)  -  (120),  the  measured  period  ratio  is  (T.  -  T„) 

o  A  D 

TBm/TAm  "  [(c  +  u)(c  "  v)/(c  "  u)(c  +  v)]1/2  *  (121) 

If  the  relativistic  relative  velocity  V  of  the  moving  points  A(u)  and  B(v)  is 
introduced,  Eq.  (121)  assumes  the  form 


TBm/TAm  “  I(c  +  V)/(c  “  (122) 

where 

V  -  (u  -  v)/(l  -  uv/c2)  .  (123) 

Since  the  measured  Doppler  frequencies  are  defined  by  v*m  -  1/T^  and  * 

1/T  ,  Eq.  (122)  formally  agrees  with  the  Doppler  effect  of  the  special  rela- 

Bm 

itivity  theory  (Whittaker,  1951;  Stratton,  1941] 
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Thus,  v*  %  v_  for  V  ^  0,  depending  on  whether  the  emitter  A  moves  towards 
Am  Bm 

(u  >  v)  or  away  from  (u  <  v)  the  moving  receiver  B.  If  A  and  B  have  the  sami 

velocity,  no  Doppler  effect  occurs,  v*  *  v_  for  u  =  v  or  V  *  0. 

Am  Bm 

It  is  remarkable  that  the  absolute  velocities  u  of  A  and  v  of  B  rela¬ 
tive  to  the  ether  frame  cannot  be  determined  by  Doppler  measurements.  The 
measured  frequency  ratio  vj^/v  gives  via  Eq.  (124)  only  the  relative  rela¬ 
tivistic  velocity  V,  since  Eq.  (123)  provides  a  manifold  of  solutions  (u,v) 
for  any  measured  value  V. 

As  another  illustration  of  the  Galilean  relativity  physics,  it  is 
recalled  that  the  Galilean  ideas  for  anisotropic  light  propagation  parallel 
(A)  and  perpendicular  (1)  to  the  ether  velocity  w  yield  for  the  fringe  shift 
A  «  At  of  the  Michelson-Morley  [1887]  experiment. 

At  -  (2/ c)  [l!1  (1  -  w2c2)-1  -  LA(1  -  v2/cV1/2]  .  (125) 

According  to  Eqs.  (Ill)  and  (112),  the  phase  difference  vanishes  exactly  since 
the  mirror  distance  is  Lorentz  contracted  whereas  the  mirror  distance  i/  is 
not  Lorentz  contracted. 

At  -  0  for  L11  -  Lq(1  -  w2/c2)1/2  ,  L1  -  Lq  .  (126) 

The  experimental  result  At  »  0  lead  Lorentz  via  Eq.  (125)  to  the  discovery  of 
the  physical  length  contraction  of  material  bodies  moving  relative  to  the  ether 
with  a  velocity  -w. 

Thus,  it  is  recognized  that  the  Galilean  concepts  of  space  and  time, 
extended  with  the  help  of  the  length  contraction  of  Lorentz  and  the  time  dila¬ 
tation  of  Builder,  explain  the  Michelson-Morley  type  interferometer  and  also 
Doppler  effect  measurements,  i.e.,  even  their  most  recent  and  highly  accurate 
versions  [Jaseja  et  al. ,  1964;  Vessot  and  Levine,  1979]. 


The  Michelson-Morley  maser  interference  experiment  of  Jaseja  et  al  [196;] 

—3  2  , 

agrees  with  Eq.  (126)  up  to  terms  of  order  10  (vq/c)  where  vq  is  the  orbital 

velocity  of  the  earth  relative  to  the  ether.  The  maser  Doppler  measurements 

of  Vessot  etal  [1972] are  in  excellent  accord  with  Eq.  (124)  since  these  permit 

to  infer  an  anisotropy  of  light  propagation  which  is  only  of  the  order  Ac/c 
-8 

<  10  .  These  investigators  believe  that  the  observed,  small  effects  are 

explainable  bymeans  of  the  general  theory  of  relativity.  The  latter  possibility 
cannot  be  discussed  here. 
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CONCLUSION 

* 

The  original  Maxwell’s  equations  for  the  ether  rest  frame  are  generalized 
to  electromagnetic  field  equations  for  arbitrary  inertial  frames,  in  which  the 
ether  is  in  a  state  of  motion,  w  4  ^3.  In  connection  with  this  theory,  we 
arrive  at  the  following  deductions. 

The  electromagnetic  field  equations  with  ether  flow  w  are  Galilei  covariant, 
and  reduce  to  the  Maxwell  equations  in  the  limit  w/c  -*■  Thus,  Maxwell's  conjec¬ 

ture  that  his  equations  hold  in  a  frame  of  reference  with  resting  ether  is  recon¬ 
firmed.  Maxwell's  equations  are  not  Galilei  covariant,  since  they  do  not  refer 
to  an  arbitrary  inertial  frame  with  moving  ether,  2^5. 

Electromagnetic  signals  propagate  isotropically  with  the  speed  of  light 
—1/2 

c  »  (pc)  relative  to  the  (moving  or  resting)  ether.  Relative  to  inertial 
frames  with  ether  flow  w,  electromagnetic  signals  propagate  anisotropically . 

The  signal  velocity  is  independent  of  the  dynamics  of  the  emitting  source,  which 

Is  typical  for  wave  propagation  in  a  carrier  (ether). 

In  inertial  frames  with  ether  flow,  the  dispersion  w  *  u(£)  of  electro¬ 
magnetic  waves  is  changed,  i.e.,  an  electromagnetic  wave  of  wave  vector  Tc  experi¬ 
ences  a  blue  -  or  red  -  shift  Aw  ■  lc*w  %  0  for  propagation  with  a  wave  vector 
component  kj|  parallel  (>)  or  antiparallel  (<)  to  the  ether  velocity  w.  In  electro¬ 
magnetic  wave  phenomena,  the  ratio  of  electric  and  magnetic  field  strengths 
|^|/|5j  *  c  +  «j|  ^  c,  is  increased  by  the  ether  flow  depending  on  whether  w^  = 

£*w/|t|  *  0. 

The  electromagnetic  ether  has  been  incorrectly  linked  with  an  "absolute 
rest  system"  since  Maxwell  [Whittaker,  1909].  The  substratum  is  a  physical  con¬ 
cept  [Ives,  1952;  Kaempffer,  1953]  whereas  absolute  rest  is  a  methaphysical  con¬ 
cept,  and  these  should,  therefore,  not  be  intermixed.  Within  the  presented 
theory,  we  have  only  assumed  that  the  ether  has  different  drift  velocities  w 


in  different  inertial  frames,  and  that  the  Galilean  relativity  principle  is 
valid  [Mansouri  and  Sexl,  1977]. 

The  electromagnetic  field  equations  (14)  -  (19)  with  substratum  flow  w  and 

their  Galilei  transformations  represent  a  field  theory  in  terms  of  absolute 

space  r  *  (x,y,z)  and  time  t  coordinates.  The  transformations  (111)  -  (112) 

indicate  how  these  are  related  to  the  measured  space  and  time  coordinates  as  a 

result  of  length  contraction  and  time  dilatation  determined  by  the  Lorentz 
■*2  2  1/2 

measure  y  =  (1  -  w  /c  )  .  As  shown,  our  theory  is  in  agreement  with  the 

Michelson-Morlev  ana  Doppler  effect  experiments. 

The  infamous  difficulties  of  the  Lorentz  covariant  and  ether-free  special 
theory  of  relativity  (applicability  to  point  particles  only,  violation  of  cau¬ 
sality  principle  for  extended  particles,  infinite  self-energy  and  self-acceler¬ 
ation  of  electrons,  infinite  zero-point  energy  density  of  vacuum, twin  paradox, 
etc)  are  removed  by  adopting  a  covariance  principle  combatible  with  an  electro¬ 
magnetic  ether.  The  abandonment  of  Lorentz  covariance  permits  the  existence  of 
extended  particles  and,  thus,  eleminates  infinite  self-energy  and  self-acceler¬ 
ation  of  the  relativistic  point  particle.  The  introduction  of  the  ether  and  a 
minimum  length  in  quantum  electrodynamics  gives  a  finite  zero-point  energy  densi 
ty  of  the  vacuum  [Winterberg,  1984],  Only  that  twin,  who  moves  relative  to  the 
ether,  experiences  an  increased  life  time,  etc. 

In  addition,  the  presented  theory  justufies  the  widespread  use  of  non- 
relativistic  electrodynamics  for  the  analysis  of  moving  conductors  in  electrical 
engineering,  magnetohydrodynamics,  and  plasma  physics  [Stratton,  1941;  Wilhelm, 
1983,1984].  We  have  considered  here  only  the  transformation  of  Maxwell's 
equations  to  inertial  frames.  The  theoretical  foundations  for  transforming 
Maxwell's  equations  to  accelerated  reference  frames  have  been  laid  by  Kottler 


[1922  a,b],  Cartan  [1924],  van  Dantzig  [1934],  Shouten  and  Hantjes  [1934],  and 
in  particular  Post  [1962,  1967,  1972,  1979].  Electromagnetic  ether  effects  in 
rotating  systems  are  important  for  Sagnac  interferometers  and  ring  lasers, 
electromagnetic  sensing  of  absolute  rotation,  and  Sagnac  type  gyroscopes  for 
navigation,  as  will  be  shown  in  further  communications. 
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Fig.  1:  Galilei  transformation  from  inertial  frame  (r\t°,w°*0)  with 

resting  substratum  to  an  inertial  frame  (r,t,w)  with  streaming 

-►  — ►  —  O  •  • 

substratum  w  *  -u,  where  moves  relative  to  with  velocity  u 
(0°  ■  0  for  t*  •  t  -  0). 


Fig.  2:  Galilei  transformation  from  inertial  frame  J(r,t,w)  to  inertial 
frame  y'(r\ t\w'),  where  j'  moves  relative  to  T  with  velocity  u 
(0  -  O'  for  t  -  t’  -  0). 

Fig.  3:  Magnetic -field  pulse  B(x,t)  (in  substratum  flow  w)  with  wave  front 
at  x  *  (W|(  +  c)t  produced  by  switching  on  a  current  j*(t)  -  JH(t) 
on  the  surface  x  -  0  of  a  conductor  o. 

Fig.  4:  Magnetic  field  pulse  B(x,t)  (in  substratum  flow  w)  with  wave  front 
at  x  »  (wj|  +  c)t  produced  by  the  motion  v(t)  ■  v0H(t)  of  a  conducting 
piston  a  transverse  to  a  magnetic  field  B0  into  the  space  x  >  0. 
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Galilei  transformation  from  Inertial  frame  £*(r\t\v#-<5)  with 
resting  substratum  to  an  Inertial  frame  £(r,t,w)  with  streaming 
substratum  w  -  -u,  where  |  moves  relative  to  with  velocity  u 
(0*  •  0  for  t*  •  t  ■  0). 


Galilei  transformation  from  inertial  frame  £(r,t,w)  to  Inertial 
frame  [(r’,t\w'),  where  ['  moves  relative  to  [  with  velocity  u 


(W|+C)t 

Fig.  3:  Magnetic  field  pulse  B(x,t)  (in  substratum  flow  w)  with  wave  front 
at  x  «  (W|j  +  c) t  produced  by  switching  on  a  current  T*(t)  *  ^H(t) 
on  the  surface  x  -  0  of  a  conductor  a. 
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Fig.  4:  Magnetic  field  pulse  B(x,t)  (in  substratum  flow  w)  with  wave  front 
at  x  •  (wjj  +  c)t  produced  by  the  motion  v(t)  -  vcH(t)  of  a  conducting 
piston  o  transverse  to  a  magnetic  field  Bc  into  the  space  x  >  0. 


